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The Causal Set approach to quantum gravity asserts that spacetime, at its smallest 
length scale, has a discrete structure. This discrete structure takes the form of a locally finite 
order relation, where the order, corresponding with the macroscopic notion of spacetime 
causality, is taken to be a fundamental aspect of nature. 

After an introduction to the Causal Set approach, this thesis considers a simple toy 
dynamics for causal sets. Numerical simulations of the model provide evidence for the 
existence of a continuum limit. While studying this toy dynamics, a picture arises of how 
the dynamics can be generalized in such a way that the theory could hope to produce more 
physically realistic causal sets. By thinking in terms of a stochastic growth process, and 
positing some fundamental principles, we are led almost uniquely to a family of dynamical 
laws (stochastic processes) parameterized by a countable sequence of coupling constants. 
This result is quite promising in that we now know how to speak of dynamics for a theory 
with discrete time. In addition, these dynamics can be expressed in terms of state models 
of Ising spins living on the relations of the causal set, which indicates how non- gravitational 
matter may arise from the theory without having to be built in at the fundamental level. 

These results are encouraging in that there exists a natural way to transform this classical 
theory, which is expressed in terms of a probability measure, to a quantum theory, expressed 
in terms of a quantum measure. A sketch as to how one might proceed in doing this is 
provided. Thus there is good reason to expect that Causal Sets are close to providing a 
background independent theory of quantum gravity. 
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Chapter 1 



Introduction to Causal Sets 



1.1 Quantum gravity in general 

The quest for a theory of quantum gravity arises from the fundamental inconsistencies of 
quantum mechanics and general relativity. Due to the philosophical differences between the 
two theories, it appears that a successful marriage of the two will occur only with a radical 
reformulation of both theories. In thinking about how to approach such a reformulation, 
one must decide which aspects of nature are fundamental, and which arise as "emergent" 
structures |57|. 



There are strong indications that nature, at its smallest length scale, has a discrete 
structure. The ultraviolet divergences of quantum field theory and the singularities of 
general relativity are two examples. The infinite "entanglement entropy" of a black hole 
in semi-classical gravity is another indication. (In fact, a continuum is not experimentally 
verifiable by any finite experiment, even in principle.) 

In the case of causal sets, two aspects are regarded as fundamental — discreteness and 
causality. The choice of causality as a fundamental notion is partly aesthetic, and partly due 
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to the great success one has in recovering other aspects of spacetime geometry from a causal 
order. The Causal Set program postulates that spacetime is a macroscopic approximation to 
an underlying discrete causal order. The other familiar properties of a spacetime manifold, 
such as its metrical geometry and Lorentzian signature, arise as "emergent" properties of 
the underlying discrete order. 

Mituo Taketani, a Japanese physicist and philosopher of science, has described the 



process of physical theory construction in terms of three distinct stages [63|. The three stages 
repeat cyclically, except after each cycle the theory is understood at a deeper level. The 
three stages, using physics terminology, may be called the phenomenological, kinematical, 
and dynamical. 

The phenomenological stage concerns itself with what physical phenomena the theory 
seeks to address. For what observational results should the theory provide explanation? 
For the Copernican model of the solar system, an example of the phenomena would be the 
retrograde motion of the planets with respect to the fixed background stars. 

In the kinematical stage one decides in what language the theory will be expressed. 
What are the basic elements of the theory, the "substance" , and how do they interrelate? 
For the case of general relativity, one chooses a manifold with a Lorentzian metric. This 
stage is very important, as in it one decides "what really exists" in nature. It determines 
the mathematical and philosophical construct on which the theory will be based. 

Once a language is set in place, and the basic constructs of the theory are chosen, then 
the remaining task is to determine how these objects behave. What are the "equations 
of motion" of the theory? In the case of General Relativity, this would be the Einstein 
equations. It selects which of the kinematical possibilities (in this example spacetime man- 
ifolds) will be realized by nature. After finishing the third stage, then one will observe new 
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phenomena which are still not explained by the theory, which then begins the phenomeno- 
logical stage of the next cycle. Each successive theory will contain those of the earlier cycles, 
generally as some limiting form of the latest theory. 

This thesis presents a step towards understanding the full quantum dynamics of space- 
time. 

Prior to this work, most knowledge about causal sets was kinematical in nature. This 
involved questions such as when a causal set is well approximated by a continuum spacetime, 
and what characteristics of a causal set one can "measure" to extract information about 
the spacetime into which it might faithfully embed. However, little was understood about 
the dynamics of causal sets. One of the primary difficulties was that most of our experience 
with dynamics was for continuum theories, where one can easily write down Lagrangians 
with differential operators. Causal sets unfortunately seemed to repel attempts to construct 
an action using simple analogy with existing continuum theories. An entirely new approach 
was needed to express the dynamics of the theory. 

There exist other discrete approaches to quantum gravity which involve a causal order- 
ing. Some examples are discussed in [28, ^, p7| ]. 



1.2 Kinematics of Causal Sets 

In the case of Causal Sets, the causal set is the kinematical "substance" of the theory. The 
kinematical stage of causal sets then concerns itself with understanding how the spacetime 
manifold emerges from the underlying discrete causal order 

Note that a causal structure is a natural choice for the "substance" of the theory, because 
it encodes all the information of a continuum spacetime (metrical geometry, topology, differ- 



^ For a more extensive introduction to Causal Sets, see ^ . 
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ential structure) save a conformal factor. Thus all that is missing is the volume information. 
However, since the theory is discrete, this arises naturally as well, from counting. So a dis- 
crete causal structure has sufficient information to encode all the kinematical framework of 
general relativity. 

1.2.1 Mathematical Definitions 

A partially ordered set, or poset, is a set S with order relation -<0 which is irreflexive 
{x -/i, X Vx € S) and transitive {x ~< y and y ~< z ^ x ^ z \/x,y,z G S). An (induced) 
subposet of a poset P is a subset P' C P whose order relation is determined by the condition 
that x -< y in P' iS X -< y va P. An interval (sometimes called an Alexandrov set or 
Alexandrov region) int[x,y] in a poset P is the induced subposet P' defined by the set of 
all elements {z\x -< z ~< y}. A causal set is defined to be a partial order for which every 
interval int[x,?/] Va;,y has finite cardinality.^ An order which obeys this condition on the 
cardinality of all intervals is called locally finite. 

In a causal set (5, -<), a pair of elements x,y £ S such that x ~< y form a relation; one 
writes that x and y are related, that x precedes y, y succeeds x {y y x), x is an ancestor 
of y, and y is a descendent or successor of x. The past of an element x £ S is the subset 
past(x) = {y £ S \ y ^ x}, i.e. the set of all ancestors of x. The past of a subset of S is 
the union of the pasts of its elements. The future of an element, future(a;), is the set of all 
descendents of x. A link is a relation for which there is no z such that x ~< z ~< y. A link is 

an irreducible relation, that is, one not implied by other relations via transitivity.^ 

^There are a number of synonyms for this mathematical object: partially ordered set, poset, (partial) 
order, transitive acyclic digraph, ... 

^At times I will employ the usual abuse of notation, referring to both the set S and the poset (5, -<) by 
the same symbol S. 

'^Causet is an abbreviated form of the term causal set. I will also use the term sub-causal set, or subcauset, 
which have the obvious meaning. 

^Links are often called "covering relations" in the mathematical literature. 
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A Hasse diagram is a pictorial representation of a poset. In the diagram, each element 
of the poset is represented as a dot. A line connects any two elements x < y related by a 
link, such that the preceding element x is drawn below the succeeding element y (as is the 
case in spacetime diagrams). 

A chain, or total order, is a set of elements for which each pair is related by -<. An 
n-chain is a chain with n elements. A path in a poset is an increasing sequence of elements, 
each related to the next by a link, i.e. it is a "chain made of links" . More precisely, a path 
C between two elements x -< y is a chain which has x as its past endpoint and y as its 
future endpoint, for which there exists no element z comparable with all elements in C such 
that X ^ z ^ y. The length of a path is its number of elements. An antichain is a set of 
elements in which no pair are related by -<. Note that this corresponds to an "achronal" 
or "acausal" set in spacetime. An n-antichain is an antichain with n elements. A maximal 
antichain A is an inextendible antichain, i.e. there does not exist any element of the poset 
which is unrelated to every element of A. A maximal antichain would correspond roughly 
to an "edgeless achronal set" , or Cauchy surface assuming global hyperbolicity. The height 
of an order is the length of the longest chain in that order (or, equivalently, the length of the 
longest path), while the width is the size of the largest antichain. An N -order is a partial 
order defined on elements. A covering graph of an order is the order's Hasse diagram, 
regarded as an undirected graph. A connected component of an order is subset of elements 
which form a connected (sub)graph in the covering graph of the order, i.e. it is a subposet 
for which any two points are connected by a set of paths. A maximal element is one which 
has no successors, i.e. an element x for which ^y £ S \ x ~< y. A minimal element is one 
which has no ancestors, i.e. an element x for which G S | y -< x. An automorphism is a 
one-to-one map of S onto itself that preserves -<. 
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1.2.2 Faithful Embedding — Correspondence with the Continuum 

Consider a (Hausdorff, paracompact) manifold M with a globally hyperbolic, time ori- 



set C = (5", -<) into a spacetime {M,gab)- We call this map a conformal embedding if 
X < y ^ f{x) £ J~{f{y)) Vx,y G S. Consider an Alexandrov region of {M,gab), 
J^{p) ^ J~{q) fo'^ every p,q £ M. The map (p is called a faithful embedding if it has 
two further properties. Firstly, the number of points n mapped into an Alexandrov region 
is equal to its spacetime volume V, up to Poisson fluctuations, i.e. the probability of getting 
n points in the region is 



where p is the density set by the fundamental length scale. Second, the mean spacing 



the characteristic length scale A over which the continuum geometry varies. (Admittedly 
these properties are not entirely well defined yet, but they provide a good heuristic notion.) 
We say that a spacetime M approximates a causal set C if there exists a faithful embedding 



The notion of faithful embedding gives a correspondence between a causal set and a 
continuum spacetime. The belief is that a given causal set embeds faithfully into a unique 
spacetime. Expressed more precisely, if (/> : C — > M and : C — > M' are two faithful 
embeddings into two spacetimes M and M' , then there exists an approximate isometry 
g : M ^ M' such that (f)' = g o (j). The isometry will only be approximate because the 

condition of faithful embedding is not sufficient to fix the entire continuum geometry — 

® Actually the requirement of a strict conformal embedding may be too stringent, as we require only a 
probabilistic reproduction of the volume information. Why should the light cones be exactly fixed? One way 
to loosen this requirement is to remove relations in the embedding with probability p. (If p <^ {R — L) / R, 
for a causet with L links and R relations, this will have negligible effect, so p must be at least ~ {R — L)/ R.) 



entable, Lorentzian, smooth metric gab- Consider a map (p : S 



M from a causal 




between points in the embedding, p 



1"^ in d dimensions, must be everywhere much less than 
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there will remain small variations in the geometry, with length scale > A, which still admit 
a faithful embedding of C. In general the precise formulation and proof of this conjecture, 
which is called the "Hauptvermutung^^ of Causal Sets, is a difficult mathematical problem. 

Given an arbitrary (past and future distinguishing) spacetime, it is easy to construct 
a causal set which will admit a faithful embedding into that spacetime. The method, 
suggestively called sprinkling, is to simply select at random, via a Poisson distribution, a 
finite set of events of the spacetime. These will correspond to the elements of the causal set. 
The causal relations of the causal set are simply those inherited directly from spacetime 
causality. 

Note that in a faithfully embedded causal set most links will then be "almost null", 
due to the Lorentz invariance of (local) spacetime. There can exist two events which are 
spatially very distant, perhaps lying in different galactic clusters, which are nevertheless 
"nearest neighbors" in the sense that they are connected by an (almost) null geodesic. 
(Also note in this connection that causal sets naturally suggest a sort of "blurring" between 
the ultraviolet and the infrared, in that a field that is highly boosted gets blueshifted, while 
at the same time it makes connection with remote parts of the universe. The extent to which 
a field can be boosted (ultraviolet cutoff) is limited by the size of the universe (infrared 
cutoff), as increasing large boosts require the existence of links connecting increasingly 
"distant" elements of the causal set.) 

So the causal set is a "discrete manifold", using the words of Riemann, from which 
the macroscopic properties of the continuum spacetime: topology, differential structure, 
metrical geometry, and causal structure, are emergent. Note that by keeping the causal 
structure as the substance of the theory, and then "counting" to get volume information, 
we have been able to recover the full metrical geometry. This succeeds because the causal 
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structure of a spacetime completely determines its metric up to a conformal factor p6| . 
The discreteness, "number" , by providing the missing volume information, restores the full 
metrical geometry. Note also that topology is restored, by inheritance from the topology 
of the manifold into which the causal set will faithfully embed. Note that the Lorentzian 
signature also arises naturally from this scheme, as the unique signature which maintains 
the distinctness of past and future directions in the causal order. 

A considerable amount of work has been done in understanding the kinematics of causal 
sets. For example it is understood how to extract spacetime dimension and proper time 
from only the discrete causal order. The dynamics, however, was much less developed before 
this work. Below is a brief sketch of some of these earlier results. 

1.2.3 Causal Set Dimension 

There are many different indicators of the dimension of a causal set, that is, ways to estimate 
the dimension d of the spacetime into which a given causal set might faithfully embed. The 
actual value of this dimension is the physical dimension. To determine directly the physical 
dimension of a causal set one would need a faithful embedding, which is difficult to achieve. 
The hope, then, is to deduce what this dimension will be by looking at certain simpler 
indicators which depend only on more easily accessible features of the causal order. 

1.2.3.1 Integral dimension indicators 

The following dimension indicators use only "conformal information" of the causal set, 
i.e. they only consider conformal embeddings in their construction, rather than faithful 
embeddings. 
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Linear dimension Linear dimension (also called combinatorial dimension) is a definition 
of dimension for a poset generally used by mathematicians. It is not quite appropriate for 
causal sets, however, as the "light cones" it uses to define the order are "square" and thus 
will not be a good dimension indicator for embeddings into Minkowski space. It is defined 
as the minimum dimension d such that the poset P can be realized as points in M" with 
the order given hy x ~< y 4^ Xi < yi Vi. In spite of its unphysical character, there are many 
results known for this notion of dimension. 



Flat conformal embedding dimension The flat conformal embedding dimension (also 
called Minkowski dimension, or causal dimension) is the minimum dimension Minkowski 
space into which the causal set can be causally embedded. Since in general one would 
expect the causet to embed into Minkowski space only locally (spacetime is locally flat), the 
more appropriate question to ask is that an (appropriately chosen) sub-causal set embed 
into Minkowski space, where the subcauset is chosen such that it is "small enough" to 
"not see the curvature" of the larger spacetime. These subsets, then, should contain the 
dimensionality information of the causet. Some conjectural bounds have been placed on 
causal dimension, in part by using results known for linear dimension. Some pixies can 
be identified, which, if present as subposets, establish a lower bound on flat conformal 



embedding dimension. See |14] and |20]. 



There are a number of difficulties with the notion of causal dimension. One is the prob- 
lem of deciding how to choose a "local" subset which may represent a "small" region of the 
spacetime. For causal sets locality is very difficult to define, because the "unit balls" of the 
Lorentzian metric contain spatially very distant points. There are indications that a notion 
of locality is maintained in causal sets [p^, p^], but it remains difficult to see how to use this 
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to construct local subsets. Another shortcoming is that it is not unreasonable to imagine 
that the dimensionality of spacetime will vary with length scale. Thus the dimensional 
information encoded in one of these subsets should have length scale dependence. If the 
appropriate subset can be chosen only after some probabilistic process of coarse graining 



(see Section 1.2.6 ), the arguments used to establish bounds on flat conformal embedding 
dimension will no longer be valid. Lastly, as alluded to in an earlier footnote, it may be too 
strict to demand an exact conformal embedding. Then, for example, one cannot use pixies 
to establish firm bounds on dimension. 

These difficulties illustrate how in general it is difficult to construct a global, exact 
indicator of dimension. The more useful notions are quasilocal, probabilistic, and fractal in 
nature, taking advantage of the volume information of a faithful embedding. 



1.2.3.2 Fractal dimension indicators 

There are other indicators of causal set dimension which depend on the volume information 
as well, i.e. they consider faithful embeddings in their definition. Consider an Alexandrov 
region mt(x,y) of a causal set C. One can derive causal set invariants by counting the 
occurrence of various substructures, such as the number of relations, number of elements, 
number of chains, number of links, etc., and compare the result with what is known for 
sprinklings into n-dimensional Minkowski space, obtaining a dimension estimate for (that 
region of) the causal set from which the region was taken. Since the number (of relations, 
say) counted will never come out exactly as that which would arise from sprinkling into an 
integral dimension Minkowski space, the dimension which obtains from this procedure will 
always be fractional, as in computing the effective dimension of fractals. 

This fractal dimension will be most meaningful if there exists a large region of C cov- 
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ered by many Alexandrov sets of (approximately) the same cardinality V, which all yield 
(approximately) the same effective dimension rics- In general different pairs of points {x, y) 
will not yield the same ncff. This will occur for a number of reasons: 

• random fluctuations. 

• sampling a different region of the causal set. 

• curvature of causet/spacetime: Because the region of spacetime enclosed within the 
interval int(x,y) in general will not be flat, the relationship between number of rela- 
tions, say, and dimension will not quite be the same as that for flat spacetime. For this 
reason it is important to choose intervals which are much smaller than the curvature 
scale. 

• scale dependence of topology, though this would vanish if we constrain V to be same 
in each region mt(x,y). 

In addition, it should be noted that the notion of "covering a region of C with many Alexan- 
drov sets" may be relying on an intuition which is not valid for causal sets. The important 
point is that in any given reference frame, almost all Alexandrov sets "look extremely null" . 
Thus it is difficult to "tile" a region with such sets, a task which looks something like 
trying to cover a two dimensional region with a collection of thickened diagonal lines. In 
general each Alexandrov set will overlap a very large number of "neighboring tiles", and 
will not cover the region in the same manner as one's intuition from Riemannian signature 
geometries suggests. 
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Volume — length scaling In d-dimensional Minkowski space, an Alexandrov set of 
"height" T (proper time between its two end points) has spacetime volume 

2V^ ^ ^vr(-^ /Ty 

d \2J {{d-l)/2)ld\2 J ^ ^ 

where is the volume of a unit d-ball. For a given Alexandrov set A = int[x, y], measure T 



by finding the longest chain connecting x and y (see sec. 1.2.4. 1| below), count the number 
of elements z in A, and invert ( |1 . 1| ) to get a dimension. This is perhaps the most obvious 
measure of dimension — to simply determine the exponent with which the volume of a 
region scales with length. 

A caveat about this scheme is the unknown coefficient nid relating length of the longest 
chain to proper time (1.3). 



Counting chains Consider an interval int[a, 6] which contains N points and C chains. 
Then in the limit of large A^, 

InA^ , , 

is a measure of the causal set's dimension. This measure is useful because it can be written 
explicitly, but is perhaps impractical because of the In In in the denominator. (( |1.2| ) follows 
from the fact that the number of chains C in an interval grows exponentially with its height 
T, while its volume (A^) grows as T"^ in d dimensions.) 

Myrheim-Meyer dimension For a causal set (5*, -<) define R to be the number of related 
pairs of elements, i.e. the number of pairs {x,y) such that x ~< y oi y ~< x Vx,y G S, and 
(following Myrheim ||4l| ) define the ordering fraction r to be the fraction of pairs of elements 
which are related, i.e. r = R/[^). A causal set which is formed by sprinkling A^ points into 
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an interval of d dimensional Minkowski space will have an ordering fraction given by 



3d!(d/2)! 
r = -- 



2 (3d/2)! 



which decreases monotonically with dimension [40|. Inverting this relationship (numerically) 
yields a fractal dimension for any given r, called the Myrheim-Meyer dimension. Since this 
measure is based on measuring a large number (R), the random fluctuations will be smaller 
than those arising from similar dimension estimators which count other quantities, making 
this a computationally efficient method to estimate causal set dimension. 

Because this measure of dimension associates a dimension to any ordering fraction r, 
it is sometimes used heuristically to specify the "dimension" of a causal set as a whole, 
without regard to whether it represents an Alexandrov set or whether the region is small 
enough not to see the spacetime curvature. 



Midpoint scaling dimension Consider an Alexandrov set int[x,y] of volume V. A 
midpoint, z, between x and y, will subdivide int[a;,y] into 3 regions int[x,z], int[z, y], and 
the remainder. If this causet is faithfully embeddable into Af^, then the volume Vr of 
the first two of these regions will be (1/2)"'^. Inverting this gives a dimension estimate of 
d = log2(V/Vfi). A convenient definition for the midpoint is to maximize the minimum of 
V{x, z) and V{z, y). 



1.2.4 Geometry 

The previous section discussed briefly how to extract some topological information from the 
discrete order. Here I mention some ideas on how to extract geometrical information. 
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1.2.4.1 Proper time 

The causal set should tell us not only whether two events are related, but "how much to 
the future" one occurred after the other. 

Consider two elements in a causal set x ~< y. The longest chain connecting them will be 
a path, which may be called a (timelike) geodesic. Note that this corresponds directly to 
the notion of timelike geodesic in continuum spacetime — it is an extremal chronological 
curve connecting x and y. For a causal set which embeds faithfully into a spacetime, 
there will usually be an extremely short path between any two related elements, e.g. one 
composed of two links, because one can always go as far out along the light cone as one 
wishes ("following a link") and likely find an element which is linked to x. Recall, however, 
that for Lorentzian geometry the appropriate (timelike) extremal path is the longest path. 
Note that in general there may be multiple longest chains passing through two elements. In 
this sense the discrete notion of geodesic departs from the continuum (in a small region), 
but still a unique path length is assigned to the pair x and y. In fact, this path length is 
proportional to the proper time interval of Minkowski spacetime, in the following sense. 

Consider a causal set which arises from sprinkling into d-dimensional Minkowski space. 
Brightwell and Bollobas ||9|, |l^] have shown that for an interval int[x,y] of volume pV , the 
length L of the longest chain satisfies L{pV)~^^'^ — > in probability, where md is an 
unknown constant which depends on the dimension of the Minkowski space. Fairly tight 
bounds can be placed on md- For d > 3 

21-1/'' 2i-l/''e™+l)l/'' , , 

and it is known that m2 = 2. Assuming that this correspondence remains in the presence 
of curvature, this provides a simple explicit method to extract timelike distances from the 
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discrete order. For simplicity, we define proper time between two related elements to be the 
length of the longest chain connecting them. 

Quite a bit is known about the fluctuations in this length as well, see e.g. [10|. Owing 



to the fact that a sprinkling into an interval of Af^ is isomorphic to a random permutation, 
and that the length of the longest increasing subsequence (which is equivalent to the length 
of the longest chain) has been studied extensively, much is known about the 2 dimensional 
I, I- 



case 



1.2.4.2 Spacelike distance 

It is difficult to construct a notion of spacelike distance on a causal set, in part because 
of the non-compactness of the Lorentz group. For example, an early proposal by 't Hooft 
p2| went essentially as follows, for M*^. For two unrelated elements {x,y), find the pair of 
elements (a, 6), such that a ^ a ~< y, b y x, and b y y, which minimizes the proper time 
(as computed in the previous section) between a and b. Unfortunately this will always turn 
out to be zero (for sprinklings of M'^ for d > 2), because there will always be a pair (a, b) 
which, by a statistical fluctuation, are linked. To see this, consider in M*^ every pair (a, b) 
where a is chosen from the intersection of the past light cones of x and y (this is where 
the maximal elements of past(x) n past(y) will lie) and b is chosen from intersection of the 
future light cones of x and y. Since these regions are noncompact, there will be an infinite 
number of (approximately) statistically independent pairs to consider, leading to certainty 
of finding a linked pair. 

However, there is a proposal for finding the spacelike distance between a maximal chain 
and a point The construction is simply this: for a geodesic 7 and a point x (the 

construction assumes that 7 is chosen such that future(a;) n 7 7^ and past(a;) n 7 7^ 0) find 
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the minimal point b in future(a;) n 7 and the maximal point a in past(x) Pi 7, and take the 
spacelike separation between 7 and x to be half the timelike distance between a and b (i.e. 
half the number of links in 7 between a and b. 

1.2.4.3 Curvature 

One way to extract curvature information from the causal set is generalize Equation ( |l.lD , 
as done in p|], to the case of non-zero curvature. From this one can extract information 



about the Ricci tensor |14|. The smaller intervals could be used to measure dimension, and 



then larger ones could be used to estimate curvature. 
1.2.5 Closed Timelike Curves 

The irreflexivity of the definition of a partial order used here is simply a convenient conven- 
tion. One could just as well have chosen to define the poset to be reflexive (x -< x Vx G 5), 
but then an added condition of acyclicity would be required: x ~< y and x ^ y ^ y x. 
Without this extra condition the order would allow cycles, corresponding (one might think) 
to closed causal curves in continuum spacetime. Note, however, that such an order would 
be sick in the sense that all the elements in such a cycle are indistinguishable from each 
other in terms of the order relation, so they might as well be regarded as a single element.[] 
In this sense causal sets "predict" that there do not exist closed causal loops in spacetime. 
Evidence indicates that the failure of closed causal curves may already be encoded into 



quantum field theory, in the form of Hawking's chronology protection conjecture |3C], which 
prohibits closed causal curves from forming via a divergence of quantum field energy density 
at a chronology horizon (a horizon which separates a region of spacetime which admits closed 



'^Perhaps this suggests that we should attach a positive integer to each element of our causal set, encoding 
the cardinality of a closed causal loop which that single element represents. 
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causal curves from one which does not). In order for our definition of a causal set to be 
consistent, something like the chronology protection conjecture must hold. 



1.2.6 Coarse graining and Scale dependent topology 

In practice it will be extremely rare that a given causal set faithfully embed into any 
spacetime. Somehow the dynamics must select four dimensional, "spacetime-like" causal 
sets. However, it is important to note that one would not expect the topology of spacetime 
to be four dimensional all the way down to the Planck scale. It is reasonable to expect 
some extra compact, Kaluza-Klein-like dimensions at small length scales. In addition, it is 
likely that even the continuum approximation itself will break down at Planck distances, 
leaving something like a "spacetime foam" . Thus some form of coarse graining will probably 
be necessary to make connection with macroscopic spacetime. However, even after coarse 
graining, it is very possible that no causal set will be precisely faithfully embeddable into a 
spacetime. A notion of an approximate embedding will likely be required, as alluded to in 

In general, there are two different approaches to coarse graining. One is to "blur" or 
"average" points. A major difficulty with this method, though, is the difficulty of maintain- 
ing Lorentz invariance, since a blurring which "looked natural" in one frame would appear 
extremely non-local in another. An alternate method is to use a decimation procedure, 
wherein some fraction of the "lattice sites" are simply ignored. This approach is easier to 
use than the blurring procedure, and it maintains Lorentz invariance. The precise method 
of coarse graining we use is simply to select some fraction of the existing elements of the 
causal set at random, ignoring the remaining elements, and inheriting the causal relations 
directly from those of the "fine-grained lattice" . The random procedure is necessary both 



18 



CHAPTER 1. INTRODUCTION TO CAUSAL SETS 



to maintain Lorentz invariance, and because of the "background- free" nature of the theory, 
which leads to an absence of any other obvious method with which to select points for 
coarse graining. Stated more precisely, a coarse grained approximation to a causet C can 
be formed by selecting a sub-causet C at random, with equal selection probability for each 
element, and with the causal order of C inherited directly from that of C, i.e. x ~< y in 
C if and only if x ~< y in C. Notice that such coarse graining is a random process, so 
from a single causet of elements, it gives us in general, not another single causet, but a 
probability distribution on the causets of m < N elements. 



For example, let us start with the 20 element causet C shown in Figure 1.1, and succes- 
sively coarse grain it down to causets of 10, 5 and 3 elements. We see that, at the largest 





or 




AGAIN: 





V 



Figure 1.1: Three successive coarse grainings of a 20-element causet 



scale shown (i.e. the smallest number of remaining elements), C has coarse-grained in this 
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instance to the 3-element "V" causet. 
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1.3 Dynamics for Causal Sets 

Because of their discrete character, many issues arise in causal set dynamics which are not 
present in the formulation of dynamics for continuum theories. Some general principles 
which need to be understood in a discrete context are general covariance, "manifoldness" 
(i.e. the emergence of a continuum at macroscopic length scales), and locality. This section 
will merely present some general issues which arise when attempting to express a dynamics 
for causal sets. A precise, detailed account follows in Chapter |3[ 

1.3.1 General covariance 

An element partial order P admits a natural labeling, which is an assignment of a non- 
negative integer 0,1,2,... A^ — 1 to each element of P such that x -< y ^ label(x) < 
label(y) Vx, y G P.^ Since a coordinate system in general relativity is simply an assignment 
of a "label" for each event of spacetime, a labeling of a causal set corresponds to a choice 
of coordinates in general relativity. The continuum analog of a natural labeling might be a 
coordinate system in which is everywhere timelike (and this in turn is almost the same 
thing as a foliation by spacelike slices). One can also consider an arbitrary labeling, which 
is an assignment of integers to the elements as above but in a manner irrespective of the 
causal order -<. This would be more closely analogous to arbitrary coordinate systems. In 
that case, there would be a well-defined gauge group — the group of permutations of the 
causet elements — and labeling invariance would signify invariance under this group, in 

analogy with diffeomorphism invariance and ordinary gauge invariance. However, we have 

*A natural labeling of an order P is equivalent to what is called a "linear extension of P" in the mathe- 
matical literature. 
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not found a useful way make use of this, and consider only natural labelings. In the context 
of causal sets, then, general covariance will translate into a statement of label independence 
of the dynamics. 

The dynamics for causal sets will be expressed as a measure defined on suitably chosen 
classes of histories, which in the context of causal sets are just the causal sets themselves. 
Generally we think of these histories as having infinite cardinality, i.e. they extend arbi- 
trarily far into the future.^ The issue of general covariance then serves to limit the sets of 
histories which have a physically meaningful measure, or equivalently, what are the physi- 
cally meaningful questions one can ask of the theory. 

In discussing this issue, it will be useful to first define the notion of a stem. A full stem 
of a causal set C is a finite subcauset S for which every element not in S succeeds a maximal 
element of S. A full stem corresponds to a completed partial history of the universe. A 
partial stem of a causet C is a finite subcauset S which contains its own past, i.e. if x G 5 
and y £ C such that y < x then y & S. 

An example of a non-generally covariant question is "What is the probability that the 
universe has a 3-chain as a full stem after N elements appear." This is not generally 
covariant because a labeling is implicit in the notion "after elements appear." A covariant 
question could ask "What is the probability that the universe has a 3-chain as a full stem, 
after the growth process runs to completion", i.e. in the limit as — > oo. It is conjectured 
that all physical questions can be expressed as a logical combination of probabilities of 
the occurrence of a given partial stem in the universe. Can a measure be formulated that 
assigns a finite answer to all such questions? 

^It is a logical possibility that the universe "ends" after some finite time, i.e. the causal set has finite 
cardinality, but we disregard this eventuality on purely metaphysical grounds. 
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1.3.2 "Manifoldness" 

Almost every causal set in no way resembles any spacetime manifold. To get a feel for 
how extreme this is, it has been shown that, in the limit of large N, the number of partial 
orders defined on elements grows as 2^^/^ (to leading order). In comparison, one may 
estimate that the number of "spacetime resembling" posets on N elements is only ~ 2^'"^^. 
Somehow the dynamics must select those that (at least at a sufficiently large length scale) 
resemble spacetimes. One might expect this selection to occur at the classical level, i.e. in 
the classical limit "dynamically preferred" causets should faithfully embed into spacetimes. 
This limit arises from the constructive interference of histories, so the "spacetime resembling 
causets" should lie at "a stationary point in the causal set action" , while those that are very 
unlike continuua should have "rapidly varying" amplitudes. Of course one needs a precise 
notion of how causets are "close to each other" to be able to speak of a stationary point. 
Our intuition comes roughly from a notion of "closeness" for Lorentzian geometries. PI 

A related question to the existence of a continuum is why the cosmological constant is 
so small. If it had it's "natural" value, of 1 in Planck units, then spacetime would have 
curvature on scales of the Planck length, meaning that there is no continuum. Thus any 
theory of quantum gravity must provide some mechanism for driving A to zero. Given that 
such a (relatively unknown) mechanism exists, causal sets provide a heuristic explanation 
for why A is not exactly zero, but fluctuates about zero with an amplitude which falls off as 
l/\/iV, where N is the volume of the universe in fundamental units |55]. Given that for the 



current era A^ ~ 10^^*^, this predicts an order of magnitude that is consistent with current 
observation. 

In addition, this problem is addressed in part by the renormalization type behavior of 



^"See ||l5|] for an interesting approach to the issue of defining a distance functional on the space of 
Lorentzian geometries. 
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the dynamics and the existence of a continuum hmit, discussed in Sections 2.3 and 2.4. 



1.3.3 Locality 

Whatever the microscopic dynamics for causal sets may look like, we expect that the con- 
tinuum approximation will be governed by an effective Lagrangian 

L,sr^-A + ^ + {R^) + ... (1.4) 

where (i?^) represents terms involving the curvature squared, etc. Dimensional consistency 
indicates that the coefficient before each term gets smaller in the expansion, since the 
curvature R has dimension 1/length^, so each R must have a coefficient of Ip, implying 
that the higher order terms have negligible contribution. However, this expression will be 
difficult to compute, even using the method of estimating curvature from counting in an 



interval mentioned in §1.2.4.3, because almost every Alexandrov set is "extremely null". 

However, the causal set is an inherently non-local object, so it is not unreasonable to 
expect that the notion of locality in the continuum will not carry over in an obvious, direct 
manner to the discrete dynamics. In fact the discrete dynamics, in its current formulation 



(see e.g. ( ^.23|) ), appears quite non-local, in that the "behavior" of a "region of the causal 
set" depends on its entire past. 



Rather than directly trying to reproduce the action of (L4) (say with some additional 
matter terms), an alternative approach is to get locality later, in the effective, continuum 
theory. Then the objective would be to use the notion of locality as a guide in choosing the 
microscopic dynamics, trying to determine what it means in this context, without worrying 
necessarily about getting an action as in ( |1.4D . If we do manage to choose an effectively 
local microscopic dynamics, then the Einstein-Hilbert action will come out "for free", given 
the dimensional arguments above (and local Lorentz invariance). For the case of causal 
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sets, this approach seems more hkely to bear fruit. 

In general it is difficult to define a discrete "lattice" which is Lorentz invariant, as 
most regular lattices that one considers are not invariant under Lorentz transformations. 
However, the set of points obtained from a random "sprinkling" into a spacetime region 
is Lorentz invariant. Taking advantage of this property of a random embedding, progress 
has been made in understanding how an effectively local action may arise in the context 
of causal sets [|2^ , |5^ . Their work shows that Lorentz invariance can be made compatible 
with locality on a lattice. 

This apparent success in combining locality, Lorentz invariance, and discreteness demon- 
strates a great advantage of causal sets. Never before have all these three aspects been 
present in a physical theory. 



Chapter 2 

Investigation of Transitive 
Percolation Dynamics 



2.1 Introduction 

The dynamics of causal sets will likely find its final expression as a quantum measure defined 
over suitably chosen classes of "histories" , where in this case a history will be simply a causal 
set.|^ One may expect, in analog with the path integral formulation of quantum mechanics, 
that the quantum measure will arise from a sum over histories, which may have a form 
similar to 



where ^ is a complex amplitude for a pair of causal sets C, C", possibly depending on a set 

of parameters {q}. A difficulty in defining the quantum measure in terms of a sum of this 

nature is that the sum would likely have to be constrained to "Schwinger histories" , which 

are pairs of histories that have the same "value" at some time "T" which is to the future 
^ Much of the text of this chapter is taken directly from [M . 





C,C' 
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of any constraint which is used to define the set of histories for which one is seeking the 
measure. Because there is no covariant notion of a time T in cosmology, and the notion 
of the "value" of a causal set at a "time T" is also difficult to define, it is difficult to see 
how to directly write down a measure of the form Instead, the quantum measure 

will probably arise via a construction analogous to that which defines the classical measure 



Even though we do not know the exact form of the summand, a question which presents 
itself is how to enumerate the causal sets which form the domain of the sum itself. This 
problem has been studied extensively, often in the context partial orders as transitive, 
acyclic, directed graphs. In particular, Kleitman and Rothschild (see also [I9|) have 



shown that, in the asymptotic limit N oo, the number of distinct orders definable on 
elements is given by 

(1 + 0{l/N))cPp^ 2^'/^+3^/2-i°g2 N/2 ^ 

where 

oo 



g 2~{j+i/2f ^ 2.1289312 



j=-oo 



for even N and 



Up — 

j=-oo 



2"^' « 2.1289368 



for odd A^. Thus, for any appreciable value of N (say N > 20), in the absence of some special 
amplitude A{C,C';{q}) which for example is zero on all but a vanishingly small fraction 



of the A^-element causets, it seems that, in practice, the sum in ( |2.1| ) must performed by a 
simulation or other approximation method. An important question then is how to sample 
the set of A^-element causets. 

There exists a simple "model" for generating partially ordered sets at random, which 
is familiar in the field of random graph theory, which we call transitive percolation. The 
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name, suggested by David Meyer [^], arises from the fact that this model can be regarded 
as a sort of one-dimensional directed percolation, where a relation i ~< j is thought of as a 
"bond" or "channel" between "sites" i and j in a one dimensional lattice (c.f. e.g. [p2| ). 
It is defined by a single real parameter p (and a non- negative integer N). To generate an 
A^-element poset at random, start with a set of elements labeled 0, 1, 2, ... — 1, and 
introduce a relation between each of the (^) pairs of elements with a probability p (with the 
element with the smaller label preceding that with the greater), where p is any real number 
in [0,1]. Since the resulting relation will not be transitive in general, form its transitive 
closure (e.g. if 2 -< 3 and 3 -< 438 then enforce that 2 -< 438). 



If transitive percolation is to be used to sample the domain of summation in ( |2.1| ), then 
we need to understand in detail the resulting distribution on the set of A^-element causal 
sets, so a weight factor can be placed into the summand to correct for the bias of the 
sampling technique. Unfortunately, it is impossible to do this, for the following reason. The 
asymptotic enumeration of A^-element orders found by Kleitman and Rothschild mentioned 
above was achieved by showing that almost all A^-orders are "3-layer" orders. (An "/-layer 
order" is one in which the set of elements is partitioned into / antichains Xi , X2 , . . . ,Xi, 
where each element of Xi precedes every element of Xj for j > i + 1, and no element of 
Xi precedes any element of Xj for i > j.) Furthermore, they found that almost all 3-layer 
orders have about A^/2 elements in X2 and about A^/4 elements in the other antichains. 
Here "almost all" means that the fraction of orders with this characteristic goes to 1 in the 
limit A^ 00. This result tells us that essentially all posets sampled will be 3-layer, so that 
the weight factor will degenerate to zero for any non-3-layer posets, which bodes ill for the 
whole approach of doing a Monte-Carlo sum. 

(In connection with generic, layered orders, Deepak Dhar |25| and Kleitman and Roth- 
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schild pq] have studied the behavior of an entropy function on these posets, S{r), where 



r is the ordering fraction defined in §1.2.3.2. They found an infinite number of first order 
phase transitions, at each of which dS/dr vanishes over a finite interval of r. The order 
parameter is the average height, which increases by one across each transition. In addition, 
they have found that, for a given r, most causets are highly time-asymmetric. The presence 
of the phase transitions suggests that there may be a continuum limit.) 

Obviously these 3-layer posets in no way resemble those which would faithfully embed 
into a spacetime. Since their number grows exponentially in A^^, one may imagine that any 
dynamics for causal sets is doomed to failure, since any Boltzmann-like weight which "only" 
grows exponentially in an extensive quantity (e.g. energy) would be insufficient to overcome 
this super-exponential entropic weight factor. Thus we have a sort of entropy catastrophe, 
forcing generic causets upon us regardless of our choice of dynamics. However, the causal 
sets generated by the transitive percolation algorithm look nothing like the generic 3-layer 
orders. If this model is to be regarded as a physical dynamics in itself, then this entropy 
catastrophe is already forestalled with this quite naive dynamical model. In fact, we can 
see that the dynamics of causal sets, being inherently non-local, would be expected to have 
an action which grows quadratically with an extensive quantity, rather than linearly. Then 
this sort of non-local action is exactly what is needed to overcome the entropic dominance 
of the generic orders. (In fact, the probability of arriving at a causal set with R related 
pairs, via the transitive percolation algorithm, grows like e^^, where /3 = — Inp acts as a 
sort of inverse temperature, c.f. [^^) Note that this situation is not so different from that 
of ordinary quantum mechanics, where the smooth paths, which form a set of measure zero 
in the space of all paths, are the ones which dominate the sum over histories in the classical 
limit. 
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One important question which has not been addressed is at what value of the Kleitman 
Rothschild result becomes valid. Enumeration of partial orders by computer shows no 
obvious tendency toward the 3-layer orders, for the meager values of N which a computer 
allows. It is possible, however unlikely, that the result will be of no consequence for causal 
sets, as it emerges only after N is much larger than will ever be needed for physically 
reasonable causal sets, say N ^ 10^^''. In any event it would be useful to have a feel for 
the "domain of validity" of this asymptotic result. 



We will see that in fact transitive percolation can be regarded as much more than just 
an algorithm to generate causal sets at random to be used in a Monte-Carlo sum over his- 
tories. It is an important special case of a generic class of "sequential growth" dynamics 
for causal sets, which will be explained in detail in Chapter ^ In particular, it has many 
appealing features, both as a model for a relatively small region of spacetime and as a cos- 
mological model for spacetime as a whole. Incidentally, it has attracted the interest of both 
mathematicians and physicists for reasons having nothing to do with quantum gravity. By 
physicists, it has been studied as a problem in the statistical mechanical field of percolation. 
By mathematicians, it has been studied extensively as a branch of random graph theory 
(a poset being the same thing as a transitive acyclic directed graph). Conversely, random 
graph theory could be construed as the theory of percolation on a complete graph. Some 
physics references on transitive percolation are [42, 24, pO, 48 1. In connection with random 



graph theory, there exist a large number of results governing the asymptotic behavior of 
posets generated in this manner |17, 13, 11, 22, 33, 53, ^. 
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2.2 Features 

2.2.1 May resemble continuum spacetime 

In computer simulations, two independent coarse- graining invariant dimension indicators, 
Myrheim-Meyer dimension and midpoint scaling dimension, tend to agree with each other, 
which is encouraging if these causal sets are to embed faithfully into spacetime with a well 
defined dimension.^ Another dimension indicator, which involves counting small subcausets 
whose frequency provides an indicator of dimension, behaves poorly. However, this measure 
of dimension is not invariant under coarse graining, so it only indicates that transitive 
percolated causal sets themselves do not directly embed faithfully into Minkowski space, 
but some appropriately chosen subcauset (e.g. coarse-grained) may still approximate a 
spacetime, which is what one would expect for a dynamics of causal sets anyway. 

In the pure percolation model, however, these dimension indicators vary with time (i.e. 
with N, as the causal set "grows") which suggests that one may wish to rescale p in such 
a way as to hold the spacetime dimension constant.^ One may ask, then, if the model can 
be generalized by having p vary with in an appropriate sense. We will see in Chapter ^ 
that something rather like this is in fact possible.^ 

2.2.2 Homogeneous 

Consider the transitive percolation algorithm described in §^, and an arbitrary element of 

a causal set generated by this model, say the one labeled 257. Its future will be some causal 
^AU these numerical calculations were performed by R. Sorkin. 

■^This is only a suggestion because these estimators neglect curvature. Transitive percolation could of 
course produce something resembling a region of a curved spacetime, such as de Sitter or anti-de Sitter. 

*It should be noted that the measure of dimension that is varying with A'^ is that of the causal set in 
its entirety, not that of a "local region". In fact, it seems that, due to the homogeneity of percolation, the 
dimension of a region will depend only on its size. Thus as N increases, the dimension associated with 
that A-element "region" (the entire causal set) changes uniformly. Then, more correctly, it is the scale 
dependence of dimension in transitive percolation which may suggest that p should vary somehow. 
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set, future(257). Because of the extreme symmetry of the transitive percolation algorithm, 
the probability distribution of future(257) will be completely independent of the structure 
of that portion of the causal set which is to the past of element 257, or spacelike to 257. 
This is clear because, regardless of what is to the past and unrelated to this element, each 
successive element will join to 257 with a fixed probability p. Thus its future will behave 
the same as that of any of the other elements. 

Therefore, the only spacetimes which a causal set generated by transitive percolation 
could hope to resemble would be (space-time) homogeneous, such as the Minkowski, de Sit- 
ter, or anti-de Sitter spacetimes. Likewise transitive percolation has no hope of resembling 
a spacetime with propagating degrees of freedom, such as gravitational waves. 

2.2.3 Time reversal invar iance 

Transitive percolation is independent of time orientation. When viewed from the perspective 
of a sequential growth dynamics, this may not be so obvious, but it is clear when viewed 
from the more static algorithm described above. 

2.2.4 Existence of a continuum limit 

Moreover, computer simulations suggest strongly that the model possesses a continuum 



limit (see Section 2.3) and exhibits scaling behavior in that limit with p scaling roughly like 
clogn/n [^]. 

2.2.5 Originary transitive percolation 

There exists another model which is very similar to transitive percolation, called "originary 
transitive percolation". It is most clearly described in terms of a "cosmological growth 
process" which is introduced in § |3.1.1| . For now, suffice it to say that the model is the same 
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as transitive percolation, except that every element (but one) must be a descendent of at 
least one other element of the causet. The net effect is that the growing causal set is required 
to have an "origin" (= unique minimum element). It turns out that originary transitive 
percolation is equivalent to ordinary transitive percolation, if one "discards" all elements 
which are not to the future of the first element. That is, if one generates a causet via 
transitive percolation with N and p, and then considers only the subcauset which contains 
the (inclusive) future of "element 0" , one obtains a model equivalent to that of originary 
transitive percolation at the same p (but of course smaller A^). 

2.2.6 Suggestive large scale cosmology 

Consider a picture of causal set cosmology which involves cosmological "bounces", where 
the causal set collapses down to a single element, and then re-expands as illustrated in 



Figure 2A. Alon et al. [|| call such an element a post, which is defined as an element which 



"post" 



Figure 2.1: Transitive percolation cosmology 

is related to every other in the causal set. In the context of percolation dynamics, they 
have proved rigorously that such cosmological bounces occur with probability 1 (if p > 0), 
from which it follows that there are infinitely many cosmological cycles, each cycle but 
the first having the dynamics of originary percolation. Then the "cosmology" of transitive 
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percolation is quite suggestive, consisting of a universe which cycles endlessly through phases 
of expansion, stasis, and contraction (via fluctuation) back down to a single element. 

Note that transitive percolation is only homogeneous on average. Thus, with this sto- 
chastic model, we see a phenomenon which does not arise in deterministic theories — a 
"locally homogeneous spacetime" which nevertheless possesses points where the universe 
contracts to "volume 1" and reexpands. Also, this means that transitive percolation cannot 
produce the entirety of an Einstein spacetime, it is only possible that its continuum limit 
yield a portion of some homogeneous spacetime. 



2.2.7 Cosmological renormalization 

In Chapter ^, the classical dynamics for causal sets will be expressed in terms of a countable 
sequence of parameters, or "coupling constants", t„. Some work by Dou pH], and more 



recently |59, describes a "cosmological renormalization" process, wherein at each cycle 
of expansion, collapse to a single point, and re-expansion, the parameters describing the 
dynamics of the causal set are "renormalized" , taking new effective values in each subsequent 
cycle. It is easy to show that the percolation dynamics are the unique "fixed points" under 
this renormalization flow, and furthermore that a large class of dynamics (choices for the 
parameters t„) converge to this fixed point as the renormalization process extends to infinity. 
Thus, if our universe is in fact described by something resembling the dynamics to be derived 
in Chapter ^, and it does undergo cycles of expansion and re-contraction, then after a long 
time it will be increasingly described by this (originary) transitive percolation dynamics. 

Furthermore, if the parameters tn take the form given in ( p.l9| ), which might be ex- 
pected for a physically reasonable dynamics, then the cosmological renormalization process 
makes an interesting "prediction" regarding the early universe. As discussed in [59|, for 
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the dynamics of ( |3.19| ), the early universe behaves as originary transitive percolation for 
a period which grows longer and longer after successive cosmological renormalizations. In 
addition, the effective p of this percolation phase diminishes as l/\/iV, with N being the 
number of elements to the past of the current cosmological epoch. Thus after a long time, 
p will be driven to an arbitrarily small value. If this is the case, then it can provide some 
explanation to the puzzle of why the universe is so large, homogeneous, and isotropic. The 
reason is that, for extremely small p, originary percolation will almost surely generate a 
"Cayley tree" , which is a tree for which each element has on average two immediate succes- 
sors. After this rapidly expanding "tree phase" the universe should make a transition into 
something resembling a spacetime which is spatially homogeneous and isotropic. 

2.2.8 Phase transitions in the early universe 

For the (non-originary) transitive percolation dynamics, it is known that there is a perco- 
lation phase transition at p ~ 1/-^^, where the causet transforms qualitatively from a large 
number of small disconnected universes^ for p < pcrit to a causet with one large universe 
and a number of much smaller disconnected universes for p > ^»crit- (This can be regarded 
as the "early universe" of transitive percolation because cosmological time can be measured 
by spacetime volume, so that small N corresponds to short or early times.) A second phase 
transition gathers the disconnected branches of the universe, leading to a single connected 
universe. This occurs near the percolation transition, at j5 ~ In N/N [^5|. In fact, this sec- 
ond phase transition occurs "at the same time" as a third, at which the fraction of elements 
to the future of element becomes very close to 1. 

There is, incidentally, still some hope of being able to reproduce the generic 3-layer 
^By the term "universe" I mean simply a connected component of the order. 
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posets by running the transitive percolation algorithm at the percolation phase transition. 
It is possible that each piece (connected component) of the causet would be a generic poset. 

2.2.9 Diffusion-like model 

The expectation value of the ordering fraction < r > of a causet generated by transitive 
percolation can be computed exactly by writing a recursion relation for the number of 
descendants of element (we'll call this element eo for short). ^ For the purposes of this 
discussion we switch to the reflexive convention for defining a partial order, i.e. replace 
the irreflexivity condition x x with reflexivity x ~< x (and re-impose acyclicity with 
x~<y~<x^x = y). Define fnik) to be the probability that, considering only those 
elements whose labels are less than n, eo has exactly k descendents. Clearly fi{k) = 6i. 
The recursion relation can be defined by noting that "at stage n" , there are two ways that 
Co can have k descendents. "At stage n — 1" either cq had k descendents, and element n 
does not "link to" any element to the future of eo, or eo had k—1 descendents, and element 
n does link to one of cq's descendents. The former event occurs with probability q^, while 
the latter with probability 1 — q^~^. Thus 

fn{k) = q'fn-lik) + (1 - g'-^)/n-l(fc " 1) • (2.2) 

The expected number of descendents of cq "at stage n" is then 

n 

E ^Uk) . 
k=l 

Because of the symmetry of the percolation algorithm, the expected number of descendents 

of element 1 at stage n is equivalent to the expected number of descendents of element at 

stage n — 1. (To see this, simply relabel the causet such that element i is relabeled i — 1.) 
® This recursion relation is due to R. Sorkin. 
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Then the expected number of relations in an A'^-element percolated causet, is J2n=i ^n, or 
switching back to the irreflexive convention, 

N 

<RiN)>= J2{xn-1). 

n=l 

This recursion relation can be evaluated quite efficiently on a computer, to yield values of 
<r>, so far for up to 2^^, accurate to numerical rounding errors. 

If this Markov process is modeled by a differential equation, the "field" / behaves as a 
wave moving at constant speed "to the right", with a diffusive character. It is possible that 
further study along these lines will lead to an understanding of the asymptotic behavior of 
this model, for example to understand the scaling behavior discussed in § ^.4[ . 

A generalization of this recursion relation is possible for the originary percolation model, 
but it is more expensive computationally than the 0{N'^) algorithm described here. 

2.2.10 Gibbsian distribution 

Transitive percolation is readily embedded in a "two-temperature" statistical mechanics 
model, and as such, happens also to be "exactly soluble" in the sense that the partition 



function can be computed exactly. Details of this model will be described in |24] 



2.3 Continuum Limit 

Here we search for evidence of a continuum limit in the transitive percolation dynamics. 
One might question whether a continuum limit is even desirable in a fundamentally discrete 
theory, but a continuum approximation in a suitable regime is certainly necessary if the 
theory is to reproduce known physics. Given this, it seems only a small step to a rigorous 
continuum limit, and conversely, the existence of such a limit would encourage the belief 
that the theory is capable of yielding continuum physics with sufficient accuracy. 
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Perhaps an analogy with kinetic theory can provide a useful illustration. In quantum 
gravity, the discreteness scale is set, presumably, by the Planck length / = (where 
K = 8ttG), whose vanishing therefore signals a continuum limit. In kinetic theory, the dis- 
creteness scales are set by the mean free path A and the mean free time r, both of which 
must go to zero for a description by partial differential equations to become exact. Corre- 
sponding to these two independent length and time scales are two "coupling constants" : the 
diffusion constant D and the speed of sound Cgound- Just as the value of the gravitational 
coupling constant Gh reflects (presumably) the magnitude of the fundamental spacetime 
discreteness scale, so the values of D and Cgound reflect the magnitudes of the microscopic 
parameters A and r according to the relations 

^ A2 A 

D ~ , Csound ~ — 

r r 

or conversely 

A , r ~ ^ . 

Csound '^sound 

In a continuum limit of kinetic theory, therefore, we must have either D ^ oi Cgound ~^ c>o. 
In the former case, we can hold Csound fixed, but we get a purely mechanical macroscopic 
world, without diffusion or viscosity. In the latter case, we can hold D fixed, but we get a 
"purely diffusive" world with mechanical forces propagating at infinite speed. In each case 
we get a well defined — but defective — continuum physics, lacking some features of the 
true, atomistic world. 

If we can trust this analogy, then something very similar must hold in quantum gravity. 
To send / to zero, we must make either G or h vanish. In the former case, we would expect 
to obtain a quantum world with the metric decoupled from non-gravitational matter; that 
is, we would expect to get a theory of quantum field theory in a purely classical background 
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spacetime solving the source-free Einstein equations. In the latter case, we would expect 
to obtain classical general relativity. Thus, there might be two distinct continuum limits of 
quantum gravity, each physically defective in its own way, but nonetheless well defined. 

For our purposes, the important point is that, although we would not expect quantum 
gravity to exist as a continuum theory, it could have limits which do, and one of these 
limits might be classical general relativity. It is thus sensible to inquire whether one of 
the classical causal set dynamics we have defined describes classical spacetimes. In the 
following, we make a beginning on this question by asking whether the special case of 
"percolated causal sets", as we will call them, admits a continuum limit at all. 

Of course, the physical content of any continuum limit we might find will depend on 
what we hold fixed in passing to the limit, and this in turn is intimately linked to how we 
choose the coarse-graining procedure that defines the effective macroscopic theory whose 
existence the continuum limit signifies. Obviously, we will want to send N ^ oo for any 
continuum limit, but it is less evident exactly how we should coarse- grain and what coarse 
grained parameters we want to hold fixed in taking the limit. Indeed, the appropriate 
choices will depend on whether the macroscopic spacetime region we have in mind is, to 
take some naturally arising examples, (i) a fixed bounded portion of Minkowski space of 
some dimension, (ii) an entire cycle of a Friedmann universe from initial expansion to final 
recollapse, or {in) an A^-dependent portion of an unbounded spacetime M that expands 
to encompass all of M as — > oo. In the sequel, we will have in mind primarily the first 
of the three examples just listed. Without attempting an definitive analysis of the coarse- 
graining question, we will simply adopt the simplest definitions that seem to us to be suited 
to this example. More specifically, we will coarse-grain by the random selection procedure 
of § 1.2.6| , and we will choose to hold fixed some convenient invariants of that sub-causal-set. 
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including the ordering fraction, which, as mentioned in §1.2.3.2, can be interpreted as the 
dimension of the spacetime region it constitutes^] As we will see, the resulting scheme has 
much in common with the kind of coarse-graining that goes into the definition of continuum 
limit in quantum field theory. For this reason, we believe it can serve also as an instructive 
"laboratory" in which this concept, and related concepts like "running coupling constant" 
and "non-trivial fixed point" , can be considered from a fresh perspective. 



2.3.1 The critical point at p = 0, N = oo 

Transitive percolation is a model of random causets which depends on two parameters, 
p G [0, 1] and G N. For a given p, the model defines a probability distribution on the set 
of A^-element causetsj^ For p = 0, the only causet with nonzero probability, obviously, is 
the A^-antichain. Now let p > 0. With a little thought, one can convince oneself that for 
A^ — > oo, the causet will look very much like a chain. Indeed it has been proved [Q] (see 
also p2|) that, as A^ — > oo with p fixed at some (arbitrarily small) positive number, r — > 1 
in probability, where r is the ordering fraction of the causal set. Note that the A^-chain has 
the greatest possible number of relations (^) , so r — > 1 gives a precise meaning to "looking 
like a chain" . 

We see that for N ^ oo, there is a change in the qualitative nature of the causet as p 
varies away from zero, and the point p = 0, A^ = oo (or p = l/N = 0) is in this sense a 
critical point of the model. It is the behavior of the model near this critical point which we 
study in detail. 

^Strictly speaking this interpretation is correct only if the causal set forms an interval or Alexandrov 
neighborhood within the spacetime, but, as mentioned earlier, the notion of Myrheim-Meyer dimension 
remains useful in this wider context. 

^Strictly speaking this distribution has gauge-invariant meaning only in the limit N ^ oo (p fixed); for 
it is only insofar as the causal set "runs to completion" that generally covariant questions can be asked. 
Notice that this limit is inherent in causal set dynamics itself, and has nothing to do with the continuum 
limit considered here, which sends p to zero as ^ oo. 
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The fact that a coarse grained causet is automatically another causet will make it easy 
for us to formulate precise notions of continuum limit, running of the coupling constant 
p, etc. In this respect, we believe that this model combines precision with novelty in such 
a manner as to furnish an instructive illustration of concepts related to renormalizability, 
independently of its application to quantum gravity. 



2.3.2 The large scale effective theory 

The transitive percolation model for generating random causal sets is a "microscopic" dy- 



namics, and the procedure described in §1.2.6 provides a precise notion of coarse graining 



(that of random selection of a sub-causal-set). On this basis, we can produce an effective 
"macroscopic" dynamics by imagining that a causet C is first percolated with N elements 
and then coarse-grained down to m < N elements. This two-step process constitutes an ef- 
fective random procedure for generating m element causets depending (in addition to m) on 
the parameters and p. In causal set theory, number of elements corresponds to spacetime 
volume, so we can interpret N/m as the factor by which the "observation scale" has been 
increased by the coarse graining. If, then, Vq is the macroscopic volume of the spacetime 
region constituted by our causet, and if we take Vq to be fixed as ^ oo, then our proce- 
dure for generating causets of m elements provides the effective dynamics at volume-scale 
Vo/m (i.e. length scale (Vo/m)^/"' for a spacetime of dimension d). 

What does it mean for our effective theory to have a continuum limit in this context? 
Our stochastic microscopic dynamics gives, for each choice of p, a probability distribution 
on the set of causal sets C with elements, and by choosing m, we determine at which 
scale to examine the corresponding effective theory. This effective theory is itself just a 
probability distribution fm on the set of m-element causets, so our dynamics will have a well 
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defined continuum limit if there exists, as ^ oo, a trajectory p = p{N) along which the 
corresponding probability distributions fm on coarse grained causets approach fixed limiting 
distributions for all m. The limiting theory in this sense is then a sequence of effective 
theories, one for each m, all fitting together consistently. (Thanks to the associative (semi- 
group) character of our coarse-graining procedure, the existence of a limiting distribution 
for any given m implies its existence for all smaller m. Thus it suffices that a limiting 
distribution fm exist for m arbitrarily large.) In general there will exist not just a single 
such trajectory p = p{N), but a one-parameter family of them (corresponding to the one 
real parameter p that characterizes the microscopic dynamics at any fixed N), and one 
may wonder whether all the trajectories will take on the same asymptotic form as they 
approach the critical point p = 1/N = 0. The asymptotic form of this trajectory has been 
studied extensively in the mathematics literature |45, 13, ^, 0, ||, ^ ^ with a variety of 



motivations, including for example the search for efficient sorting algorithms. 

Consider first the simplest nontrivial case, m = 2. Since there are only two causal 
sets of size two, the 2-chain and the 2-antichain, the distribution /2 that gives the "large 
scale physics" in this case is described by a single number which we can take to be /2(l), 
the probability of obtaining a 2-chain rather than a 2-antichain. (The other probability, 
/2(»» ), is of course not independent, since classical probabilities must add up to unity.) 
Interestingly enough, the number /2(l) has a direct physical interpretation in terms of the 
Myrheim- Meyer dimension of the fine-grained causet C. Indeed, it can be seen that /2(l) 
is nothing but the expectation value of what we called above the ordering fraction r of 
C (an argument explaining why this is so follows in the next section). But the ordering 
fraction, in turn, determines the Myrheim-Meyer dimension d that indicates the dimension 
of the Minkowski spacetime M'^ (if any) in which C would embed faithfully as an interval 
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pO| , |4l| ]. Thus, by coarse graining down to two elements, we are effectively measuring a 
certain kind of spacetime dimensionality of C. In practice, we would not expect C to embed 
faithfully without some degree of coarse-graining, but the original r would still provide a 
good dimension estimate since it is, on average, coarse-graining invariant. 

As we begin to consider coarse-graining to sizes m > 2, the degree of complication grows 
rapidly, simply because the number of partial orders defined on m elements grows rapidly 
with m. For m = 3 there are five possible causal sets: I, V , I , , and Thus 
the effective dynamics at this "scale" is given by five probabilities (so four free parameters). 
For m = 4 there are sixteen probabilities, for m = 5 there are sixty-three, and for m = 6, 7 
and 8, the number of probabilities is respectively 318, 2045, and 16999. 



2.3.3 Evidence from simulations 

In order that a continuum limit exist, it must be possible to choose a trajectory for p 
as a function of N so that the resulting coarse-grained probability distributions, /i, /2, 
/s, . . . , have well defined limits as — > oo. To study this question numerically, one can 



simulate transitive percolation using the algorithm described in Section |11|, while choosing 
p so as to hold constant (say) the m = 2 distribution /2 (/i being trivial). Because of 
the way transitive percolation is defined, it is intuitively obvious that p can be chosen to 
achieve this, and that in doing so, one leaves p with no further freedom. (Observe that 
< r >= /2( I ) is when p = 0, 1 when p = 1, and increases smoothly and monotonically 
with p. Thus for any choice of < r >G [0, 1] there must a p which yields that < r >, and since 
f2{'* ) = 1-/2(1)) the entire distribution /2.) The decisive question then is whether, along 
the trajectory thereby defined, the higher distribution functions, /s, fi, etc. all approach 
nontrivial limits. 
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As mentioned above, holding /2 fixed is the same thing as holding fixed the expectation 
value <r > of ordering fraction r = R/{^). To see in more detail why this is so, consider 
the coarse-graining that takes us from the original causet On oi N elements to a causet C2 
of two elements. Since coarse- graining is just random selection, the probability /2(l) that 
C2 turns out to be a 2-chain is just the probability that two elements of Cat selected at 
random form a 2-chain rather than a 2-antichain. In other words, it is just the probability 
that two elements of Cat selected at random are causally related. Plainly, this is the same 
as the fraction of pairs of elements of Cat such that the two members of the pair form a 
relation x y or y ~< x. Therefore, the ordering fraction r equals the probability of getting 
a 2-chain when coarse graining Cn down to two elements; and /2(l) =<t>, as claimed. 

This reasoning illustrates, in fact, how one can in principle determine any one of the 
distributions fm by answering the question, "What is the probability of getting this partic- 
ular m-element causet from this particular A^-element causet if you coarse grain down to m 
elements?" To compute the answer to such a question starting with any given causet Cn, 
one examines every possible combination of m elements, counts the number of times that 
the combination forms the particular causet being looked for, and divides the total by (^). 
The ensemble mean of the resulting abundance, as we will refer to it, is then fm{C)^ where 
^ is the causet for which one is looking. In practice, of course, we normally use a more 
efficient counting algorithm than simply examining individually all (^) subsets of Cjq ■ 

2.3.3.1 Histograms of 2-chain and 4-chain abundances 

As explained in the previous subsection, the main computational problem, once the random 
causet has been generated, is determining the number of subcausets of different sizes and 
types. To get a feel for how some of the resulting "abundances" are distributed, we start 
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Figure 2.2: Distribution of number of relations for N = 4096, p = 0.01155 



6,822,375 



by presenting a couple of histograms. Figure 2.2 shows the number R of relations obtained 



from a simulation in which 15,260 causal sets were generated by transitive percolation with 
p = 0.01155, = 4096. Visually, the distribution is Gaussian, in agreement with the fact 
that its "kurtosis" 



{x — x)'^ / {x — x)^ 



of 2.993 is very nearly equal to its Gaussian value of 3 (the over-bar denotes sample mean). 
In these simulations, p was chosen so that the number of 3-chains was equal on average to 
half the total number possible, i.e. the "abundance of 3-chains", (number of 3-chains) / (^) , 
was equal to 1/2 on average. The picture is qualitatively identical if one counts 4-chains 
rather than 2-chains, as exhibited in Fig. p.3|. 



(One may wonder whether it was to be expected that these distributions would appear 
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Figure 2.3: Distribution of number of 4-chains for N = 4096, p = 0.01155 

to be so normal. If the variable in question, here the number of 2-chains R or the number 
of 4-chains (C4, say), can be expressed as a sum of independent random variables, then the 
central limit theorem provides an explanation. So consider the variables Xij which are 1 if 
i ~< j and zero otherwise. Then R is easily expressed as a sum of these variables: 

R = ^ Xij 

i<j 

However, the Xij are not independent, due to transitivity. Apparently, this dependence is 
not large enough to interfere much with the normality of their sum. The number of 4-chains 
C4 can be expressed in a similar manner 

XijX jjf^Xkl . 

i<j<k<l 

and similar remarks apply.) 

Let us mention that for values of p sufficiently close to or 1, these distributions will 
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appear skew. This occurs simply because the numbers under consideration (e.g. the number 
of ?n-chains) are bounded between zero and (^) and must deviate from normahty if their 
mean gets too close to a boundary relative to the size of their standard deviation. Whenever 
we draw an error bar in the following, we will ignore any deviation from normality in the 
corresponding distribution. 

Notice incidentally that the total number of 4-chains possible is ("^^^ ^) = 11, 710, 951, 848, 960. 
Consequently, the mean 4-chain abundance^ in our simulation is only ^ "^jiq ^g^{ ^gjg ^g^Q = 

f\ "TOO '7QO 

0.234, a considerably smaller value than the mean 2-chain abundance of <r>= 'mo96\ = 



0.802. This was to be expected, considering that the 2-chain is one of only two possible 
causets of its size, while the 4-chain is one out 16 possibilities. (Notice also that 4-chains 
are necessarily less probable than 2-chains, because every coarse-graining of a 4-chain is a 
2-chain, whereas the 2-chain can come from every 4-element causet save the 4-antichain.) 

2.3.3.2 Trajectories of p versus N 

The question we are exploring is whether there exist, for oo, trajectories p = p{N) 

along which the mean abundances of all finite causets tend to definite limits. To seek such 
trajectories numerically, we will select some finite "reference causet" and determine, for a 
range of N, those values of p which maintain its mean abundance at some target value. If a 
continuum limit does exist, then it should not matter in the end which causet we select as 
our reference, since any other choice (together with a matching choice of target abundance) 
should produce the same trajectory asymptotically. We would also anticipate that all the 
trajectories would behave similarly for large A^, and that, in particular, either all would 

lead to continuum limits or all would not. In principle it could happen that only a certain 

^Occasionally I will write simply "abundance", in place of "mean abundance", assuming the average is 
obvious from context. 
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subset led to continuum limits, but we know of no reason to expect such an eventuality. 
In the simulations reported here, I have chosen as our reference causets the 2-, 3- and 
5-chains. I have computed six trajectories, holding the mean 2-chain abundance fixed at 
1/2, 1/3, and 1/10, the mean 3-chain abundance fixed at 1/2 and .0814837, and the mean 
5-chain abundance fixed at 1/2. For N, I have used as large a range as available computing 
resources allowed. 
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Figure 2.4: Ordering fractions as a function of p for = 2048 



Before discussing the trajectories as such, let us have a look at how the mean 2-chain 
abundance < r > (i.e. the mean ordering fraction) varies with p for a fixed of 2048, as 
exhibited in Figure (Vertical error bars are displayed in the figure but are so small 
that they just look like horizontal lines. The plotted points were obtained from the exact 
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calculation for the ensemble average <r> discussed in §2.2.£, so the errors come only from 



floating point roundoff. The fitting function used in Figure 2.4 will be discussed in |48|. 
As one can see, < r > starts at for p = 0, rises rapidly to near 1 and then asymptotes 
to 1 at p = 1 (not shown). Of course, it was evident a priori that < r > would increase 
monotonically from to 1 as p varied between these same two values, but it is perhaps 
noteworthy that its graph betrays no sign of discontinuity or non-analyticity (no sign of a 
"phase transition"). To this extent, it strengthens the expectation that the trajectories we 
find will all share the same qualitative behavior as — > oo. 
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Figure 2.5: Flow of the "coupling constant" p as — > oo (six trajectories) 



The six trajectories we have simulated are depicted in Fig. |2.5| . P^ A higher abundance 
of m-chains for fixed m leads to a trajectory with higher p. Also note that, as observed 



^"Notice that the error bars are shown rotated in the legend. This will be the case for all subsequent 
legends as well. 
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above, the longer chains require larger values of p to attain the same mean abundance, 
hence a choice of mean abundance =1/2 corresponds in each case to a different trajectory. 
The trajectories with <r> held to lower values are "higher dimensional" in the sense that 
<r >= 1/2 corresponds to a Myrheim-Meyer dimension of 2, while <r >= 1/10 corresponds 
to a Myrheim-Meyer dimension of 4. Observe that the plots give the impression of becoming 
straight lines with a common slope at large N. This tends to corroborate the expectation 
that they will exhibit some form of scaling with a common exponent, a behavior reminiscent 
of that found with continuum limits in many other contexts. This is further suggested by 
the fact that two distinct trajectories (/2(l) = 1/2 and /3(i) = .0814837), obtained by 
holding different abundances fixed, seem to converge for large N. 
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Figure 2.6: Six trajectories approaching the critical point at p = 0, N = oo 



By taking the abscissa to be 1/A^ rather than log2A^, we can bring the critical point 
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to the origin, as in Fig. 2.6. The lines which pass through the data points there are just 
sphnes drawn to aid the eye in fohowing the trajectories. Note that the curves tend to 
asymptote to the p-axis, suggesting that p fahs off more slowly than This suggestion 

is corroborated by more detailed analysis of the scaling behavior of these trajectories, as 
will be discussed in [Ell. 



2.3.3.3 Flow of the coarse-grained theory along a trajectory 

We come finally to a direct test of whether the coarse-grained theory converges to a limit 
as iV ^ oo. Independently of scaling or any other indicator, this is by definition the 
criterion for a continuum limit to exist. I have examined this question by means of sim- 
ulations conducted for five of the six trajectories mentioned above. In each simulation I 
proceeded as follows. For each chosen N, I experimentally found a p sufficiently close to 
the desired trajectory. Having determined p, I then generated a large number of causets 
by the percolation algorithm described in Section |TT|. (The number generated varied from 
64 to 40,000.) For each such random causet, I computed the abundances of the different 
m-element (sub)causets under consideration (2-chain, 3-chain, 3-antichain, etc), and com- 
bined the results to obtain the mean abundances we have plotted here, together with their 
standard errors. (The errors shown do not include any contribution from the slight inac- 
curacy in the value of p used. Except for the 3- and 5-chain trajectories these errors are 
negligibly small.) 

To compute the abundances of the 2-, 3-, and 4-orders for a given causet, I randomly 
sampled its four-element subcausets, counting the number of times each of the sixteen 
possible 4-orders arose, and dividing each of these counts by the number of samples taken 
to get the corresponding abundance. As an aid in identifying to which 4-order a sampled 
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subcauset belonged I used the following invariant, which distinguishes all of the sixteen 
4-orders, save two pairs. 

/(5) = n(2 + |past(x)|) 

Here, past(3;) = {y £ S\y -< x} is the exclusive past of the element x and |past(x)| is its 
cardinality. Thus, we associate to each element of the causet, a number which is two more 
than the cardinality of its exclusive past, and we form the product of these numbers (four, 
in this case) to get our invariant. (For example, this invariant is 90 for the "diamond" 
poset, ^.) 

The number of samples taken from an N element causet was chosen to be y^2(^), on the 
grounds that the probability to get the same four element subset twice becomes appreciable 
with more than this number of samples. Numerical tests confirmed that this rule of thumb 
tends to minimize the sampling error, as seen in Figure |2.7| . F^ 

Once one has the mean abundances of all the 4-orders, the mean abundances of the 
smaller causets can be found by further coarse graining. By explicitly carrying out this 
coarse graining, one easily deduces the following relationships: 

/3(i) = /4(i) + ^(/4(r)+/4(k))+^/4(i.) + ^(/4(I^) + /4(h))+^/4(I>) 
/3(V) = ^/4(r) + ^/4(l^) + ^/4(I>) + ^/4(Y) + ^/4(Y.) + ^/4(H) + ^/4(M) 
/3(I.) = ^/4(i.) + ^(/4(I^)+/4(h))+^(/4(Y.)+/4(A.))+/4(ll) + ^/4(l..) 

/3(A) = ^/4(k) + ^/4(h) + ^/4(I>) + ^/4(A) + ^/4(A.) + ^/4(H) + ^/4(M) 

/3(—) = ^(/4(Y)+/4(A))+^(/4(Y.)+/4(A.))+^/4(l..)+/4(....) 

^^The errors depicted in Fig. were found by generating 100 causets by transitive percolation, and for 
each one performing the indicated number of samples of 4 element subcausets (with replacement), counting 
the fraction of times that the diamond arose. The errors reported are the (square root of the) variance of 
the mean of this quantity over the 100 causets. 
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Figure 2.7: Reduction of error in estimated diamond abundance with increasing number of 
samples 

/2(I) = /3(i) + ^(/3(V)+/3(A))+^/3(I.) 
/2(..) = 1-/2(1) 



In the first six equations, the coefficient before each term on the right is the fraction of 
coarse- grainings of that causet which yield the causet on the left. 

In Figures U, U, and |l0|, we exhibit how the coarse-grained probabilities of all 
possible 2, 3, and 4 element causets vary as we follow the trajectory along which the 
coarse-grained 2-chain probability /2(l) = ^ is held at 1/2. By design, the coarse-grained 



probability for the 2-chain remains flat at 50%, so Figure 2.8 simply shows the accuracy 



with which this was achieved. (Observe the scale on the vertical axis.) Notice that, since 
72(1) and /2(*» ) must sum to 1, their error bars are necessarily equal. (The standard 
deviation in the abundances decreases with increasing N. The "blip" around log2 = 9 
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Figure 2.8: Flow of the coarse-grained probabilities fm for m = 2. The 2-chain probability 
is held at 1/2. 



occurs simply because I generated fewer causets at that and larger values of N to reduce 
computational costs.) 

The crucial question is whether the probabilities for the three and four element causets 
tend to definite limits as tends to infinity. Several features of the diagrams indicate that 
this is indeed occurring. Most obviously, all the curves, except possibly a couple in Figure 



2.1C, appear to be leveling off at large N. But we can bolster this conclusion by observing 



in which direction the curves are moving, and considering their interrelationships. 



For the moment let us focus our attention on Figure 2.9. A priori there are five coarse- 
grained probabilities to be followed. That they must add up to unity reduces the degrees 
of freedom to four. This is reduced further to three by the observation that, due to the 
time-reversal symmetry of the percolation dynamics, we must have /3( V) = f3{j\), as 
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Figure 2.9: Flow of the coarse-grained probabilities fm for m = 3. The 2-chain probability 
is held at 1/2. 



duly manifested in their graphs. Moreover, all five of the curves appear to be monotonic, 
with the curves for j\, \^ and ••• rising, and the curves for \ and I , falling. If we 
accept this indication of monotonicity from the diagram, then first of all, every probability 
/3(^) must converge to some limiting value, because monotonic bounded functions always 
do; and some of these limits must be nonzero, because the probabilities must add up to 1. 
Indeed, since /3( V) and f3{/\) are rising, they must converge to some nonzero value, 
and this value must lie below 1/2 in order that the total probability not exceed unity. In 
consequence, the rising curve /3(»»») must also converge to a nontrivial probability (one 
which is neither nor 1). Taken all in all, then, it looks very much like the m = 3 coarse- 
grained theory has a nontrivial N ^ oo limit, with at least three out of its five probabilities 
converging to nontrivial values. 
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Figure 2.10: Flow of the coarse-grained probabilities fm for m = 4. The 2-chain probability 
is held at 1/2. 



Although the "rearrangement" of the coarse-grained probabilities appears much more 



dramatic in Figure 2.10, similar arguments can be made. Excepting initial "transients", it 



seems reasonable to conclude from the data that monotonicity will be maintained. From 
this, it would follow that the probabilities for ^ and (which must be equal by time- 
reversal symmetry) and the other rising probabilities, and all approach 
nontrivial limits. The coarse-graining to 4 elements, therefore, would also admit a contin- 
uum limit with a minimum of 4 out of the 11 independent probabilities being nontrivial. 

To the extent that the m = 2 and m = 3 cases are indicative, then, it is reasonable 
to conclude that percolation dynamics admits a continuum limit which is non-trivial at all 
"scales" m. 
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Figure 2.11: Flow of the coarse-grained probabilities fm for m = 2. The 3-chain probability 
is held at 0.0814837. 



The question suggests itself, whether the flow of the coarse-grained probabilities would 
differ qualitatively if we held fixed some mean abundance other than that of the 2-chain. In 



Figures 2.11 , 2.12 , and 2.15 , we display results obtained by fixing the 3-chain abundance (its 
value having been chosen to make the abundance of 2-chains be 1/2 when N = 2^^). Notice 



in Figure 2.11 that the abundance of 2-chains varies considerably along this trajectory, 



whilst that of the 3-chain (in Figure 2.12) of course remains constant. Once again, the 
figures suggest strongly that the trajectory is approaching a continuum limit with nontrivial 
values for the coarse-grained probabilities of at least the 3-chain, the "V" and the "A" (and 
in consequence, of the 2-chain and 2-antichain) . 

All the trajectories discussed so far produce causets with an ordering fraction r close to 
1/2 for large N. As mentioned earlier, r = 1/2 corresponds to a Myrheim-Meyer dimension 
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Figure 2.12: Flow of the coarse-grained probabilities fm for m = 3. The 3-chain probability 
is held at 0.0814837. 



of two. Figures |2.14| and |2.15| show the results of a simulation along the "four dimensional" 
trajectory defined by r = 1/10. (The value r = 1/10 corresponds to a Myrheim- Meyer 
dimension of 4.) Here the appearance of the flow is much less elaborate, with the curves 
arrayed simply in order of increasing ordering fraction, ••• and •••• being at the top and 
i and (imperceptibly) | at the bottom. As before, all the curves are monotone as far as 
can be seen. Aside from the intrinsic interest of the case d = 4, these results indicate that 
our conclusions drawn for d near 2 will hold good for all larger d as well. 



Figure 2.16 displays the flow of the coarse-grained probabilities from a simulation in 
the opposite situation where the ordering fraction is much greater than 1/2 (the Myrheim- 
Meyer dimension is down near 1.) Shown are the results of coarse-graining to three element 
causets along the trajectory which holds the 3-chain probability to 1/2. Also shown is the 
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Figure 2.13: Flow of the coarse-grained probabilities fm for m = 4. The 3-chain probabihty 
is held at 0.0814837. 



2-chain probability. The behavior is similar to that of Figure 2.14 , except that here the 
coarse-grained probability rises with the ordering fraction instead of falling. This occurs 
because constraining /3(I) to be 1/2 generates rather chain- like causets whose Myrheim- 
Meyer dimension is in the neighborhood of 1.34, as follows from the approximate limiting 
value /2(l) ~ 0.8. The slow, monotonic, variation of the probabilities at large N, along 
with the appearance of convergence to non-zero values in each case, suggests the presence 
of a nontrivial continuum limit for r near unity as well. 



Figures p. 17 and 2.18 present the results of a final set of simulations, the only ones 
we have carried out which examined the abundances of causets containing more than four 
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Figure 2.14: Flow of the coarse-grained probabilities fm for m = 3. The 2-chain probability 
is held at 1/10. 

elements. In these simulations, the mean 5-chain abundance /5(5-chain) was held at 1/2, 
producing causets that were even more chain-like than before (Myrheim-Meyer dimension 



~ 1.1). In Figure 2.17, we track the resulting abundances of all A;-chains for k between 2 and 
7, inclusive. (We limited ourselves to chains, because their abundances are relatively easy to 



determine computationally.) As in Figure 2. If, all the coarse-grained probabilities appear 
to be tending monotonically to limits at large A^. In fact, they look amazingly constant 
over the whole range of A^, from 5 to 2^^. One may also observe that the coarse-grained 
probability of a chain decreases markedly (and almost linearly over the range examined!) 
with its length, as one might expect. It appears also that the A;-chain curves for k ^ 5 are 
"expanding away" from the 5-chain curve, but only very slightly. Figure 2.18| displays the 
flow of the probabilities for coarse-grainings to four elements. It is qualitatively similar to 



Figures 2.14-2.16, with very flat probability curves, and here with a strong preference for 
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Figure 2.15: Flow of the coarse-grained probabilities fm for m = 4. The 2-chain probability 
is held at 1/10. Only those curves lying high enough to be seen distinctly have been labeled. 



causets having many relations over those having few. Comparing Figures 2.18 and 2.15 with 



Figures 2.13 and 2.10, we observe that trajectories which generate causets that are rather 
chain-like or antichain-like seem to produce distributions that converge more rapidly than 
those along which the ordering fraction takes values close to 1/2. 

In the way of further simulations, it would be extremely interesting to look for contin- 



uum limits of some of the more general dynamical laws discussed in § p.3.5| . In doing so, 
however, one would no longer have available (as one does have for transitive percolation) a 
very fast (yet easily coded) algorithm that generates causets randomly in accord with the 
underlying dynamical law. Since the sequential growth dynamics of Chapter ^ is produced 
by a stochastic process defined recursively on the causal set, it is easily mimicked algorithmi- 
cally; but the most obvious algorithms that do so are too slow to generate efficiently causets 
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Figure 2.16: Flow of the coarse-grained probabilities fm for m = 3. The 3-chain probability 
is held at 1/2. 



of the size we have discussed in this paper. Hence, one would either have to devise better 
algorithms for generating causets "one off", or one would have to use an entirely different 
method to obtain the mean abundances, like Monte Carlo simulation of the random causet. 



2.3.4 Concluding Comments 

Transitive percolation is a discrete dynamical theory characterized by a single parameter p 



lying between and 1. Regarded as a stochastic process as described in ^3.1.1, it describes 
the steady growth of a causal set by the continual birth or "accretion" of new elements. If 
we limit ourselves to that portion of the causet comprising the elements born between step 
A'^o and step A'^i of the stochastic process, we obtain a model of random posets containing 
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N = Ni- No elements. 

Because the underlying process is homogeneous, this model does not depend on A'o or A'^i 
separately, but only on their difference. It is therefore characterized by just two parameters 
p and A^. One should be aware that this truncation to a finite model is not consistent 
with discrete general covariance (c.f. § 3.2.2 and § 1.3.1| ), because it is the subset of elements 
with certain labels that has been selected out of the larger causet, rather than a subset 
characterized by any directly physical condition. Thus, we have introduced an "element of 
gauge" and we hope that we are justified in having neglected it. That is, we hope that the 
random causets produced by the model we have actually studied are representative of the 
type of suborder that one would obtain by percolating a much larger (eventually infinite) 
causet and then using a label- invariant criterion to select a subset of N elements. 
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Figure 2.18: Flow of the coarse-grained probabilities fm for m = 4. The 5-chain probability 
is held at 1/2. 

Leaving this question aside, let us imagine that our model represents an interval (say) in 
a causet C underlying some macroscopic spacetime manifold. With this image in mind, it 
is natural to interpret a continuum limit as one in which N ^ oo while the coarse-grained 
features of the interval in question remain constant. We have made this notion precise by 
defining coarse-graining as random selection of a suborder whose cardinality m measures 
the "coarseness" of our approximation. A continuum limit then is defined to be one in 
which tends to oo such that, for each finite m, the induced probability distribution fm 
on the set of m-element posets converges to a definite limit, the physical meaning being 
that the dynamics at the corresponding length-scale is well defined. Now, how could our 
model fail to admit such a limit? 



In a field-theoretic setting, failure of a continuum limit to exist typically means that the 
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coarse-grained theory loses parameters as the cutoff length goes to zero. For example, A^^ 
scalar field theory in 4 dimensions depends on two parameters, the mass fi and the coupling 
constant A. In the continuum limit, A is lost (i.e. it necessarily vanishes, because it is driven 
to infinity at finite cutoff if one tries to hold its renormalized value fixed), although one 
can arrange for /i to survive. (At least this is what most workers believe occurs.) Strictly 
speaking, one should not say that a continuum limit fails to exist altogether, but only that 
the limiting theory is poorer in coupling constants than it was before the limit was taken. 
Now in our case, we have only one parameter to start with, and we have seen that it does 
survive as — > oo since we can, for example, choose freely the m = 2 coarse-grained 
probability distribution f2- Hence, we need not fear such a loss of parameters in our case. 

What about the opposite possibility? Could the coarse-grained theory gain parameters 
in the N ^ oo limit, as might occur if the distributions were sensitive to the fine details 
of the trajectory along which N and p approached the "critical point" p = 0, N = oo?^ Our 
simulations showed no sign of such sensitivity, although we did not look for it specifically. 
(Compare, for example. Figure 2^ with Figure |1| and with |1|.) 



A third way the continuum limit could fail might perhaps be viewed as an extreme form 
of the second. It might happen that, no matter how one chose the trajectory p = p{N), some 
of the coarse-grained probabilities fmiO oscillated indefinitely as A^ — > oo, without ever 
settling down to fixed values. Our simulations leave little room for this kind of breakdown, 
since they manifest the exact opposite kind of behavior, namely monotone variation of all 
the coarse-grained probabilities we "measured". 

Finally, a continuum limit could exist in the technical sense, but it still could be ef- 
fectively trivial (once again reminiscent of the Xc/)^ case — if you care to regard a free 



^^For example, at a water-ice phase transition, the state is specified by the pressure, temperature, and 
density, while away from the transition the pressure and temperature suffice. Thus the theory picks up a 
parameter in the limit of approaching the phase transition. 
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field theory as trivial). Here triviality would mean that all — or almost all — of the 
coarse-grained probabilities fmiO converged either to or to 1. Plainly, we can avoid 
this for at least some of the fm{0- -^^^ example, we could choose an m and hold either 
/m("i-chain) or /m(?72-antichain) fixed at any desired value. (As p — > 1, fm{m-cham) — > 1 
and /m(m-antichain) ^0; as p — > 0, the opposite occurs.) However, in principle, it could 
still happen that all the other fm besides these two went to in the limit. (Clearly, they 
could not go to 1, the other trivial value.) Once again, our simulations show the opposite 
behavior. For example, we saw that /3( V) increased monotonically along the trajectory 



of Figure 2.9 



Moreover, even without reference to the simulations, we can make this hypothetical 
"chain-antichain degeneracy" appear very implausible by considering a "typical" causet C 
generated by percolation for ^ 1 with p on the trajectory that, for some chosen m, holds 
/m(?Ti-chain) fixed at a value a strictly between and 1. Then our degeneracy would insist 
that /m(?7T--antichain) = 1 — a and fm{x) = for all other x- I^ut this would mean that, in 
a manner of speaking, "every" coarse-graining of C to m elements would be either a chain 
or an antichain. In particular the causet I , could not occur as a subcauset of C; whence, 
since I , is a subcauset of every m-element causet except the chain and the antichain, C 
itself would have to be either an antichain or a chain. But it is absurd that percolation 
for any parameter value p other than and 1 would produce a "bimodal" distribution such 
that C would have to be either a chain or an antichain, but nothing in between. It seems 
likely that similar arguments could be devised against the possibility of similar, but slightly 
less trivial trivial continuum limits, for example a limit in which fmix) would vanish unless 
X were a disjoint union of chains and antichains. 

Putting all this together, we have persuasive evidence that the percolation model does 
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admit a continuum limit, with the hmiting model being nontrivial and described by a single 
"renormalized" parameter or "coupling constant". Furthermore, the associated scaling 
behavior one might think to find in such a case is also present, as will be discussed further 



m 



ii 



However, the question remains as to whether this continuum resembles a spacetime 
manifold, or is something more pathological. Do the causal sets yielded by the percolation 
dynamics resemble genuine spacetimes? Based on the meager evidence available at the 



present time, for example that mentioned in §2.2.1 and §2.2.2 we can only answer "it is 
possible" . 

Finally, there is the ubiquitous issue of "fine tuning" or "large numbers" . In any contin- 
uum situation, a large number is being manifested (an actual infinity in the case of a true 
continuum) and one may wonder where it came from. In our case, the large numbers were 
and A^. For N, there is no mystery: unless the causal set "ends" at some point, N is 
guaranteed to be as large as desired. But why should p be so small? Here we can appeal 



to the preliminary results of Dou mentioned in § 2.2.7 , which state that under cosmological 



renormalization, certain physically reasonable dynamics of causal sets are driven toward an 

— 1/2 

effective percolation-like phase with a value of p that scales like Nq , where Nq is the 
number of elements of the causet preceding the most recent "bounce" . Since this is sure 
to be an enormous number if one waits long enough, p is sure to become arbitrarily small 
if sufficiently many cycles occur. The reason for the near flatness of spacetime — or if 
you like for the large diameter of the contemporary universe — would then be just that 
the underlying causal set is very old — old enough to have accumulated, let us say, 10^'^'^ 
elements in earlier cycles of expansion, contraction and re-expansion. 
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2.4 Scaling 



The trajectories of Figure 2.5 can be fit numerically to explore the scaling behavior of this 
model. Doing so seems to suggest an asymptotic functional form of log N/N, which agrees 
with the prediction of |45]. Not only does this evidence for scaling support the conclusion 



of the existence of a continuum limit, but the asymptotic form can also be used to suggest 
a choice for the parameters tn of the classical dynamics (of ( ^.21[) , say) which may produce 
spacetimes with a dimension which is constant over many length scales. Details will appear 
in [|8|. 



Chapter 3 



Classical Dynamics of Sequential 
Growth 

3.1 Introduction 

As illustrated in Chapter |l], the causal set approach to quantum gravity has experienced 
considerable progress in its kinematic aspects.^ For example, there exist natural extensions 
of the concepts of proper time and spacetime dimensionality, which take us a significant 
way toward an answer to the question, "When does a causal set resemble a Lorentzian 
manifold?". The dynamics of causal sets (the "equations of motion"), however, has not 
been very developed to date. One of the primary difficulties in formulating a dynamics for 
causal sets is the sparseness of the fundamental mathematical structure. When all one has 
to work with is a discrete set and a partial order, even the notion of what we should mean 
by a dynamics is not obvious. 

Traditionally, one prescribes a dynamical law by specifying a Hamiltonian to be the 
^Much of the text of this chapter is taken directly from ]49[. 
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generator of the time evolution. This practice presupposes the existence of a continuous 
time variable, which we do not have in the case of causal sets. Thus, one must conceive 
of dynamics in a more general sense. Considerable progress can be made by envisaging 
evolution as a process of stochastic growth to be described in terms of the probabilities (in 
the classical case, or more generally the quantum measures in the quantum case [pq]) of 



forming designated causal sets. As mentioned in § |1.3.1| , the dynamical law will be a rule 
which assigns probabilities to suitable classes of causal sets (a causal set being the "history" 
of the theory in the sense of "sum-over-histories" ) . One can then use this rule — technically 
a probability measure — to ask physically meaningful questions of the theory. For example 
one could ask "What is the probability that the universe possesses the diamond poset as a 
partial stem?". 

Why are we interested in a classical dynamics for causal sets, when our ultimate aim 
is a quantum theory of gravity? One obvious reason is that the classical case, being much 
simpler, can help us to get used to a relatively unfamiliar type of dynamical formulation, 
bringing out the pertinent physical issues and guiding us toward physically suitable condi- 
tions to place on the theory. Is there, for example, an appropriate form of causality that 
we can impose? Should we attempt to express the theory directly in terms of gauge in- 
variant (labeling independent) quantities, or should we follow precedent by enforcing gauge 
invariance only at the end? Some of these issues are well illustrated with the theories we 
construct herein. 

One of the best reasons to be interested in a classical dynamics for causal sets is that 
quantum gravity must possess general relativity as a classical limit. Thus general relativity 
should be described as some type of effective classical dynamics for causal sets, and one 
may hope that the relevant dynamical law will be among the family delineated here. (One 
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can't be certain this will occur, because general relativity, as a continuum theory, seems 
most likely to arise as an effective theory for coarse-grained causal sets, rather than directly 
as a limit of the microscopic discrete theory, and there is no guarantee that this effective 
theory will have the same form as the underlying exact one.) 

A question commonly asked of the causal set program is "How could nongravitational 
matter arise from only a partial order?" . One obvious answer is that matter can emerge as a 
higher level construct via the Kaluza-Klein mechanism [p9| ] , but this possibility has nothing 
to do with causal sets as such. The theory developed in this chapter suggests a different 
mechanism, in that is possible to rewrite the theory in such a way that the dynamics appears 
to arise from a kind of "effective action" for a field of Ising spins living on the relations of 
the causal set. A form of "Ising matter" is thus implicit in what would seem at first sight 
to be a purely "source-free" theory. 

3.1.1 Sequential growth 

The dynamics which we will derive can be regarded as a process of "cosmological accretion" 
or "growth" . At each step of this process an element of the causal set comes into being as 
the "offspring" of a definite set of the existing elements - the elements that form its past. 
The phenomenological passage of time is taken to be a manifestation of this continuing 
growth of the causet. Thus, we do not think of the process as happening "in time" but 
rather as "constituting time" , which means in a practical sense that there is no meaningful 
order of birth of the elements other than that implied by the relation -< . 

In order to define the dynamics, however, we will treat the births as if they happened in 
a definite order with respect to some fictitious "external time". In this way, we introduce 
an element of "gauge" into the description of the growth process which we will have to 
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compensate by imposing appropriate conditions of "gauge invariance" . This fictitious order 
of birth can be represented as a natural labeUng of the elements, that is, a labeling by 
integers 0, 1, 2, 3, . . . which are compatible with the causal order (recall definition of natural 
labeling in § 1.3.1| ). The relevant notion of gauge invariance (which we will call "discrete 



general covariance") is then captured by the statement that the labels carry no physical 
meaning. We discuss this more extensively later on. 

It is helpful to visualize the growth of the causal set in terms of paths in a poset V 
of finite causal sets. (Thus viewed, the growth process will be a sort of Markov process 
taking place in V.) Each finite causet (or rather each isomorphism equivalence class of 
finite causets) is one element of this poset. If a causet can be formed by accreting a single 
element to a second causet, then the former (the "child") follows the latter (the "parent") 
in V and the relation between them is a link. Drawing V as a Hasse diagram of Hasse 
diagrams, we get Fig. |3.1| . (Of course this is only a portion of the infinite diagram; it 
includes all the causal sets of fewer than five elements and 8 of the 63 five element causets. 



The "decorations" on some of the transitions in Fig. 3A are for later use.) Any natural 
labeling of a causet C £ V determines uniquely a path in V beginning at the empty causet 
and ending at C. Conversely, any choice of upward path through this diagram determines a 
naturally labeled causet, or rather a set of them, since inequivalent labelings can sometimes 
give rise to the same path in We want the physics to be independent of labeling, so 
different paths in V leading to the same causet should be regarded as representing the same 
(partial) universe, the distinction between them being "pure gauge" . 

The causal sets which can be formed by adjoining a single maximal element to a given 
causet will be called collectively a family. The causet from which they come is their parent, 



^We could restore uniqueness by "resolving" each link Ci ^ C2 oi V into the set of distinct embeddings 
i : Ci —> C2 that it represents. Here, two embeddings count as distinct iflt no automorphism of the child 
relates them (cf. the discussion of the Markov sum rule below). 
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i q =1 

Figure 3.1: The poset of finite causets 

and they are siblings of each other. Each one is a child of the parent. The child formed 
by adjoining an element which is to the future of every element of the parent will be called 
the timid child. The child formed by adjoining an element which is spacelike to every other 
element will be called the gregarious child. A child which is not the timid child will be 
called a bold child. (At times I may not be careful to distinguish between a child and the 
transition probability leading to that child, assuming that the intended meaning is obvious 
from context.) 

Each parent-child relationship in V describes a 'transition' C — > C, from one causal 
set to another induced by the birth of a new element. The past of the new element (a 
subset of C) will be referred to as the precursor set of the transition (or sometimes just the 
"precursor of the transition"). Normally, this precursor set is uniquely determined up to 
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automorphism of the parent by the (isomorphism equivalence class of the) child, but (rather 
remarkably) this is not always the case.|^ The symbol C„ will denote the set of causets with 




Figure 3.2: Inequivalent precursors lead to same child 

n elements, and the set of all transitions from C„ to Cn+i will be called stage n. 

As just remarked, each parent-child transition corresponds to a choice of partial stem 
in the parent (the precursor of the transition). Since there is a one-to-one correspondence 
between partial stems and antichains, a choice of child also corresponds to a choice of 
(possibly empty) antichain in the parent, the antichain in question being the set of maximal 
elements of the precursor. Note also that the new element will be linked to each element of 
this antichain. 



3.1.2 Some examples 



To help clarify the terminology introduced in the previous section, we give some examples. 



The 20 element causet of Figure ^ was generated by the stochastic dynamics described 
herein, with the choice of parameters given by Equation ( |3.19D below (with t = 1). In the 
copy of this causet on the left, the past of element a is highlighted. Notice that since we 



■^An example of the latter situation is shown in Figure 3.2, in which a parent causal set C (whose elements 
are the 10 squares) can undergo a transition to a child causal set C' by adjoining a new element, which is 
shown as one of the two circles. Observe that a new element at either circle will lead to the same 11 element 
causet, but there is no automorphism of the parent which maps one precursor into the other. 
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Past of element 'a' 



A partial stem 



Figure 3.3: An example of a ('typical'?) 20 element causal set 

use the irreflexive convention for the order, a is not included in its own past. In the the 
copy on the right, a partial stem of the causet is highlighted. 

Figure 3.4 shows , and its children. The timid child is Of, and the gregarious child is 



Cc- The precursor set leading to the transition to Cd is shown in the ellipse. An example of 
an automorphism of Ca is the map a <^ c,b <^ d (the other elements remaining unchanged) . 

3.2 Physical requirements on the dynamics 



The dynamics of transitive percolation, which was introduced in § |ll| , can be expressed as 
a growth dynamics of the sort presented in the previous section, by stating that each new 
element forges a causal bond independently with each existing element with probability 
p S [0,1]. (Any causal relation implied by transitivity must then be added in as well.) 
However, this is only one special case drawn from a much larger universe of possibilities. 
As preparation for describing these more general possible dynamical rules, let us consider 
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timid cliild 



child 



greganous 
child 




precursor 
set 



parent 



Figure 3.4: A family 



the growth-sequence of a causal set universe. 

First element '0' appears (say with probability one, since the universe exists). Then 
element '1' appears, either related to '0' or not. Then element '2' appears, either related to 
'0' or '1', or both, or neither. Of course if 1 ;^ and 2 y 1 then 2 by transitivity. Then 
element '3' appears with some consistent set of ancestors, and so on and so forth. Because 
of transitivity, each new element ends up with a partial stem of the previous causet as its 
precursor set. The result of this process, obviously, is a naturally labeled causet (finite if we 
stop at some finite stage, or infinite if we do not) whose labels record the order of succession 



of the individual births. For illustration, consider the path in fig. 3.1 delineated by the 



heavy arrows. Along this path, element '0' appears initially, then element '1' appears to 
the future of element '0', then element '2' appears to the future of element '0', but not to 
the future of '1', then element '3' appears unrelated to any existing element, then element 
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'4' appears to the future of elements '0', '1' (say, or '2', it doesn't matter) and '3', then 
element '5' appears (not shown in the diagram), etc. 

Let us emphasize once more that the labels 0, 1, 2, etc. are not supposed to be physically 
significant. Rather, the "external time" that they record is just a way to conceptualize the 
process, and any two birth sequences related to each other by a permutation of their labels 
are to be regarded as physically identical. 

So far, we have been describing the kinematics of sequential growth. In order to define 
a dynamics for it, we may give, for each n-element causet C, the transition probability 
from it to each of its possible children. Equivalently, we give a transition probability for 
each partial stem within C. We wish to construct a general theory for these transition 
probabilities by subjecting them to certain natural conditions. In other words, we want to 
construct the most general (classically stochastic) "sequential growth dynamics" for causal 
sets.0 In stating the following conditions, we will employ the terminology introduced above. 

3.2.1 The condition of internal temporality 

By this imposing sounding phrase, we mean simply that each element is born either to the 
future of, or unrelated to, all existing elements; that is, no element can arise to the past of 
an existing element. 

We have already assumed this tacitly in describing what we mean by a sequential growth 

dynamics. An equivalent formulation is that the labeling induced by the order of birth must 

be natural, as defined above. The logic behind the requirement of internal temporality is 
*By choosing to specify our stochastic process in terms of transition probabilities, we have assumed in 
effect that the process is Markovian. Although this might seem to entail a loss of generality, the loss is 
only apparent, because the condition of discrete general covariance introduced below would have forced the 
Markov assumption on us, even if we had not already adopted it. 
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that all physical time is that of the intrinsic order defining the causal set itself. For an 
element to be born to the past of another would be contradictory: it would mean that an 
event occurred "before" another which intrinsically preceded it. 

3.2.2 The condition of discrete general covariance 

As we have been emphasizing, the "external time" in which the causal set grows (equiv- 
alently the induced labeling of the resulting poset) is not meant to carry any physical 
information. We interpret this in the present context as being the condition that the net 
probability of forming any particular n-element causet C is independent of the order of birth 
we attribute to its elements. Thus, if 7 is any path through the poset V of finite causal sets 
that originates at the empty causet and terminates at C, then the product of the transition 
probabilities along the links of 7 must be the same as for any other path arriving at C. 
(So general covariance in this setting is a type of path independence). We should recall 
here, however, that, as observed earlier, a link in V can sometimes represent more than one 
possible transition. Thus our statement of path-independence, to be technically correct, 
should say that the answer is the same no matter which transition (partial stem) we select 
to represent the link. Obviously, this immediately entails that all such representatives share 
the same transition probability. 

We might with justice have required here conditions that are apparently much stronger, 
including the condition that any two paths through V with the same initial and final end- 
points have the same product of transition probabilities. However, it is easy to see that 
this already follows from the condition stated]^ We therefore do not make it part of our 

definition of discrete general covariance, although we will be using it crucially. 

^If 7 does not start with the empty causet Co, but at Cs, we can extend it to start at Co by choosing any 
fixed path from Co to Cs- Then different paths from Cs to the end-point Ce correspond to different paths 
between Co and Ce , and the equality of net probabilities for the latter implies the same thing for the former. 
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Finally, it is well to remark here that just because the "arrival probability at C" is 
independent of path/labeling, that does not necessarily mean that it carries an invariant 
meaning. On the contrary a statement like "when the causet had 8 elements it was a chain" 
is itself meaningless before a certain birth order is chosen. This, also, is an aspect of the 
gauge problem, but not one that functions as a constraint on the transition probabilities 
that define our dynamics. Rather it limits the physically meaningful questions that we can 
ask of the dynamics. Technically, we expect that our dynamics (like any stochastic process) 
can be interpreted as a probability measure on a certain cr-algebra, and the requirement 
of general covariance will, in addition providing a constraint on the transition probabilities 
of the growth process, serve to select the subalgebra of sets whose measures have direct 
physical meaning. 

3.2.3 The Bell causality condition 

The condition of "internal temporality" may be viewed as a very weak type of causality 
condition. The further causality condition we introduce now is a discrete analog of the 
statement that no influence can propagate faster than light. This condition is quite strong, 
being similar to that from which one derives Bell's inequalities. We believe that such a 
condition is appropriate for a classical theory, and we expect that some analog will be 
valid in the quantum well. (On the other hand, we would have to abandon Bell 

causality if our aim were to reproduce quantum effects from a classical stochastic dynamics, 
as is sometimes advocated in the context of "hidden variable theories" . Given the inherent 
non-locality of causal sets, there is no logical reason why such an attempt would have to 
fail.) 

The physical idea behind our condition is that events occurring in some part of a causal 
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set C should be influenced only by the portion of C lying to their past. In this way, the 
order relation constituting C will be causal in the dynamical sense, and not only in name. 
In terms of our sequential growth dynamics, we make this precise as the requirement that 
the ratio of the transition probabilities leading to two possible children of a given causet 
depend only on the triad consisting of the two corresponding precursor sets and their union. 

Thus, let C — > Ci designate a transition from C € C„ to Ci G Cn+i, and similarly for 
C — > C2. Then, the Bell causality condition can be expressed as the equality of two ratios]^ 

probjC ^ Ci) ^ probjB ^ ^1) 

prob{C C2) prob{B ^ B2) ^ ' ' 

where B G Cm-, m < n, is the union of the precursor set of C ^ Ci with the precursor 
set of C — > C2, Si G Cm+i is B with an element added in the same manner as in the 
transition C ^ Ci, and B2 G Cm+i is B with an element added in the same manner as in 
the transition C ^ C2.[| (Notice that if the union of the precursor sets is the entire parent 
causet, then the Bell causality condition reduces to a trivial identity.) 

To clarify the relationships among the causets involved, it may help to characterize the 
latter in yet another way. Let ei be the element born in the transition C — > Ci and let 62 
be the element born in the transition C C2. Then Cj = C U {e,/} {i = 1, 2), and we have 
B = (past ei) U (past 62) and Bi = B U {ej} {i = 1,2). 

By its definition. Bell causality relates ratios of transition probabilities belonging to 
one "stage" of the growth process to ratios of transition probabilities belonging to previous 
stages. For illustration, consider the case depicted in fig. |3.5| . The precursor Pi of the 
transition C ^ Ci contains only the earliest (minimum) element of C, shown in the figure 
as a pattern- filled dot. The precursor P2 of C ^ C2 contains as well the next earliest 



In writing (3.1), we assume for simplicity that both numerators and both denominators are nonzero. 

'^Recall that the precursor set of the transition C — » Ci is the subposet of C that lies to the past of the 
new element that forms Ci. 
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Figure 3.5: Illustrating Bell causality 



element, shown as a (different pattern)-filled dot. The union of the two precursors is thus 
B = PiU P2 = P2- The elements of C depicted as open dots belong to neither precursor. 
Such elements will be called spectators. Bell causality says that the spectators can be 
deleted without affecting relative probabilities. Thus the ratio of the transition probabilities 



of Figure 3.5 is equal to that of Figure 3.6 



,0 
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Figure 3.6: Illustrating Bell causality - spectators do not affect relative probability 



3.2.4 The Markov sum rule 

As with any Markov process, we must require that the sum of the full set of transition 
probabilities issuing from a given causet be unity. However, the set we have to sum over 
depends in a subtle manner on the extent to which we regard causal set elements as "dis- 
tinguishable" . Heretofore we have identified distinct transitions with distinct precursor sets 
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of the parent. In doing so, we have in effect been treating causet elements as distinguish- 
able (by not identifying with each other, precursor sets related by automorphisms of the 
parent), and this is what we shall continue to do. Indeed, this is the counting of children 
used implicitly by transitive percolation, so we keep it here for consistency. With respect to 



the diagram of Figure 3.1, this method of counting has the effect of introducing coefficients 
into the sum rule, equal to the number of partial stems of the parent which could be the 
precursor set of the transition. For the transitions depicted there, these coefficients (when 
not one) are shown next to the corresponding arrow .0 

These sum-rule coefficients admit an alternative description in terms of embeddings of 
the parent into the child (as a partial stem). Instead of saying "the number of partial stems 
of the parent which could be the past of the new element" , we could say "the number of 
order preserving injective maps from the parent onto partial stems of the child, divided by 
the number of automorphisms of the child". To see this, let e be such an embedding of 
parent P into child C. 

P ^C 

Given that the child has precisely one more element than the parent, the injective map e 
singles out an element as the unique member of C not belonging to the image of e. This 
element can be regarded as the new element that arises from the transition P ^ C . The 
past of this element (in P), n(e), will then be a possible precursor of this transition 

n(e) = e-^(past(C\e(P))) . (3.2) 

However, the maps e overcount the number of precursors by the number of automorphisms 

®One might describe the result of setting these coefficients to unity as the case of "indistinguishable 
causet elements" . A preliminary investigation suggests that in this case a dynamics with a richer structure 
obtains: instead of the transition probability depending only on the size of the precursor set and the number 
of its maximal elements, it is sensitive to more details of the precursor set's structure. However, since the 
dynamics of transitive percolation does not satisfy this modified sum rule, it appears difficult to derive a 
closed form expression for the transition probabilities in this case. 
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of the child, since each e composed with an automorphism a G Aut(C) yields a new map 
f = ao e which corresponds to the same precursor. To prove this, we must show that 

U{e)=U{f)^f = aoe 

To prove insert ( p^ ) for / 

n(/) = (Qoe)-i(past(C\a(e(P)))) 

= {ae)-\j>ast{a{C\e{P))))f\ 

= (ae)-^Q(past(C\e(P))) 

= e-^a-^a(past(C\e(P))) 

= n(e). 

The third equality holds because a is an automorphism. To prove write / as /3 o e, for 
some map (3 : C ^ C. We then have 

e-i(past(C7\e(P))) = /-i(past(C \ /(P))) 

= e-V(past(C\/3(e(P)))) 
= e-V-'(past(/?(C\e(P)))). 

For this last equality to hold in general, (3 must be an automorphism. 

\ oi{X) = oi{C \ X) as sets, for any set X C C. An element 

y ea(C\X) 

^ a-^y eC\X 
<^ a^^y <f. X 
*^yi ol{X) , 

which is exactly the condition that y £ C \ a(X). 
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3.3 The general form of the transition probabiUties 

We seek to derive a general prescription which gives, consistent with our requirements, the 
transition probability from an element of Cn to an element of Cn+i- To avoid having to deal 
with special cases, we will assume throughout that no transition probability vanishes. Thus 
the solution we find may be termed "generic" , but not absolutely general. 

3.3.1 Counting the free parameters 

A theory of the sort we are seeking provides a probability for each transition, so without 
further restriction, it would contain a free parameter for every possible antichain of every 
possible (finite) causet. We will see, however, that the requirements described above in 



Section 3.2 drastically limit this freedom. 



Lemma 1 There is at most one free parameter per family. 

Proof: Consider a parent and its children (the set of possible transitions from the parent). 
A Bell causality equation will relate any pair of transitions whose union of precursors is 
not the entire parent causet, since the remaining elements will then be spectators, whose 
removal provides such a relation. Since the precursor for the gregarious child is empty, the 
complement of its union with any other partial stem (save that of the timid child, whose pre- 
cursor is the entire parent) will be non-empty, resulting in a Bell causality equation relating 
the pair of transition probabilities. Thus every child, except the timid child, participates 
in a Bell causality equation with the gregarious child. (See also the proof of Lemma |5| in 
Appendix ^.) Since Bell causality equates ratios, all these transitions are determined up 
to an overall factor. This leaves two free parameters for the family. The Markov sum rule 
gives another equation, which exhausts itself in determining the probability of the timid 
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child. Hence precisely one free parameter per family remains after Bell causality and the 
sum rule are imposed. □ 

Lemma 2 The probability to add a completely disconnected element (the "gregarious child 
transition") depends only on the cardinality of the parent causal set. 

Proof: Consider an arbitrary causet A, with a maximal element e, as indicated in Figure 




Figure 3.7: Equality of "gregarious child" transitions 



3.7. Adjoining a disconnected element to A produces the causet B. Then, removing e 
from B leads to the causet C, which can be looked upon as the gregarious child of the 
causet D = A\{e}. Adding another disconnected element to C leads to a causet E with (at 
least) two completely disconnected elements. The lower case letter attached to each arrow 
represents the corresponding transition probability. Now, by general covariance, 



ax = bw 



and by Bell causality, 
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(the disconnected element in C acts as the spectator here). Thus 



ax = bw = ay =^ x = y 

(RecaU that we have assumed that no transition probabihty vanishes.) Repeating our 
deductions with C in the place of A in the above argument (and another maximal element 
/ in the place of e), we see that y = z, where z is the probability for the transition from F 
to G (which has at least three completely disconnected elements) as shown. Continuing in 
this way until we reach the antichain An shows finally that x = qn, where we define qn as 
the transition probability from the n-antichain to the (n + l)-antichain. Since our starting 
causet A was not chosen specially, this completes the proof. □ 

If our causal sets are regarded as entire universes, then a gregarious child transition 
corresponds to the spawning of a new, completely disconnected universe (which is not to 
say that this new universe will not connect up with the existing universe in the future). 
Lemma |2| proves that the probability for this to occur does not depend on the internal 
structure of the existing universe, but only on its size, which seems eminently reasonable. 
In the sequel, we will call this probability g„. 

With Lemmas [l| and ^, we have reduced the number of free parameters (since every 
family has a gregarious child) to 1 per stage, or what is the same thing, to one per causal 
set element. In the next sections we will see that no further reduction is possible based 
on our stated conditions. Thus, the transition probabilities qn can be identified as the free 
parameters or "coupling constants" of the theory. They are, however, restricted further by 
inequalities that we will derive below. 
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3.3.2 The general transition probability in closed form 

Given the the remaining transition probabilities (for the non- gregarious children) are 
determined by Bell causality and the sum rule, as we have seen. Here we derive an expression 
in closed form for an arbitrary transition probability in terms of causet invariants and the 
parameters g„. 

3.3.2.1 Mathematical form of transition probabilities 

We use the following notation: 
a„ an arbitrary transition probability from Cn to Cn+i 

Pn a transition whose precursor set is not the entire parent ('bold' transition) 

7„ a transition whose precursor set is the entire parent ('timid' transition) 

Notice that the subscript n here refers only to the number of elements of the parent causet; 
it does not exhibit which particular transition of stage n is intended. A more complete 
notation might provide a, (3 and 7 with further indices to specify both the parent causet 
and the precursor set within the parent. 
We also set go = 1 by convention. 

Lemma 3 Each transition probability an of stage n has the form 



where the are integers depending on the individual transition in question. 

Proof: This is easily seen to be true for stage 1. Assume it is true for stage n — 1. Consider 
a non-timid transition probability (3n of stage n. Bell causality gives 




(3.3) 



•n 



(3.4) 



qn-i 
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where an-i is an appropriate transition probability from stage n — 1. (Bell causality has 
the property that not all spectators have to be removed in writing the right hand side of 
the equation. Any partial stem of the subcauset of spectators may be "kept" in selecting 
which transition probabilities to place on the right hand side.) So by induction 

n— 1 n— 1 

For a timid transition probability 7„, we use the Markov sum rule: 

7n = 1 - E (3-^) 

j 

where j labels the possible bold transitions (i.e. the set of proper partial stems of the 
parent).!^ But then, substituting (3.5) yields immediately 



n-1 ^ n-1 e.. 

7n = i-EE^^n = i-E^' 

j i=0 i=0 



which we clearly can put into the form (3.3) by taking = — J2j ^ji fo'^ i < n and ^„ = 1. 

□ 

3.3.2.2 Another look at transitive percolation 

The transitive percolation model introduced in Chapter is consistent with the four con- 



ditions of Section |3.2| . To see this, consider an arbitrary causal set C„ of size n. Recall 
that, expressed in terms of a sequential growth process, transitive percolation states that 
at each stage the new element joins to each pre-existing element with probability p, with 
extra relations added to insure transitivity. Then the transition probability a„ from Cn to 
a specified causet Cn+i of size n + 1 is given by 

a„=p-(l-p)"— (3.7) 



'Of course, more than one stem will in general correspond to the same link in V. If we redefined j to run 



over links in V, then (15.6) would read 7„ = 1 — ^ , XjPnj, where Xj is the "multiplicity" of the jth link. 
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where m is the number of maximal elements in the precursor set and -oj is the size of the 
entire precursor set. (This becomes clear if one recalls how the precursor set of a newborn 
element is generated in transitive percolation: first a set of ancestors is selected at random, 
and then the ancestors implied by transitivity are added. From this, it follows immediately 
that a given stem S* C C„ results from the procedure iff (i) every maximal element of S 
is selected in the first step, and (ii) no element of Cn\S is selected in the first step.) In 
particular, we see that the "gregarious transition" will occur with probability g„ = (7", 
where q = 1 — p. 

Now consider our four conditions. Internal temporality was built in from the outset, as 
we know. Discrete general covariance is seen to hold upon writing the net probability of a 
given Cn explicitly in terms of causet invariants (writing it in "manifestly covariant form" , 
c.f. the discussion in § |3.5| ) as 

Pr(C„) = I^(C„)p^g(2)-^ 

where L is the number of links in Cn, R the number of relations, and W the number of 
(natural) labelings of C„. (To see how this arises, note that each of the links in the causal 
set have to be "put in by hand", i.e. a link will never enter during the transitive closure 
stage of the algorithm, so the L links occur with probability p^. Furthermore, each of the 
(2) — R non-relations of the causet must have not been "selected" during pre-transitive 
closure stage of the algorithm, each of which occurs with probability q = 1—p. Finally, the 
transitive percolation algorithm generates labeled posets, so any given unlabeled causet Cn 
can arise in W{Cn) different ways.) 

To see that transitive percolation obeys Bell causality, consider an arbitrary parent 
causet. The transition probability to a given child is exhibited in eq. ( |3.7| ). Consider two 
different children, one with (m, ro) = (mi,tZ7i) and the other with {m, -07) = (1712,^2). Bell 
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causality requires that the ratio of their transition probabihties be the same as if the parent 
were reduced to the union of the precursor sets of the two transitions, i.e. it requires 

pm2qn—-W2 pm2qn'—m2 

where n' is the cardinahty of the union of the precursor sets of the two transitions. Thus, 
Bell causality is satisfied by inspection. 

Finally, the Markov sum rule is essentially trivial. At each stage of the growth process, a 
preliminary choice of ancestors is made by a well-defined probabilistic procedure, and each 
such choice is mapped uniquely onto a choice of partial stem. Thus the induced probabilities 
of the partial stems sum automatically to unity. 

3.3.2.3 The general transition probability 

In the previous section we have shown that transitive percolation produces transition prob- 
abilities ( |3.7| ) consistent with all our conditions. By equating the right hand side of ( ^.T] ) 
to the general form ( |3.3D of Lemma ^, we can solve for the and thus obtain the general 
solution of our conditions: 

n I 

an = Y.^^-qn= " p)""^ = (1 " (Z)™?""^ 

Expanding the factor (1 — g)'", and using the fact that g„ = for transitive percolation, 
we get 

f. = (-r-f '".). 

\w — ij 

So an arbitrary transition probability in the general dynamics is, according to (|3.3|) 
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Noting that the binomial coefficients are zero for w — i ^ {0..m}, and rearranging the 
indices, we obtain 



This form for the transition probabihty exhibits its causal nature particularly clearly: except 
for the overall normalization factor g^, depends only on invariants of the associated 
precursor set. 

3.3.3 Inequalities 

Since the are classical probabilities, each must lie between and 1, and this in turn 
restricts the possible values of the Here we show that it suffices to impose only one 
inequality per stage; all the others (two per child) then follow. More precisely, what we 
show is that, if g„ > for all n, and if a„, > for the "timid" transition from the n-antichain, 
then all the he in [0, 1]. This we establish in the following two "Claims". 

Claim In order that all the transition probabilities an fall between and 1, it suffices that 
each timid transition probability 6e > 0. 

Proof: As described in the proofs of lemmas |^ and ^, each bold transition (of stage n) is 
given (via Bell causality) by 



where m is some natural number less than n. The g's are positive. So if the probabilities 
of the previous stages are positive, then the bold probabilities of stage n are also positive. 
It follows by induction that all but the timid transition probabilities are positive (since 
CKo = 50 = 1 obviously is). But for the timid transition of each family, we have 






(3.9) 
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where each /3j is positive. If any of the (3i is greater than one, jn will obviously be negative. 
Also ( |3.9| ) plainly cannot be greater than one. Consequently, if we require that 7„ be 
positive, then all transition probabilities in the family will be in [0, 1]. □ 

In a timid transition, the entire parent is the precursor set, so w = n. The inequalities 
constraining each probability of a given family to be in [0, 1] therefore reduce to the sole 
condition 

Claim The most restrictive inequality of stage n is the one arising from the n-antichain, 
i.e. the one for which m = n. All other inequalities of stage n follow from this inequality 
and the inequalities for smaller n. 



Proof: Assume that we have, for m = n, 



fc=0 



n\ 1 



klq, 



> 0. 



n—k 



Add to this the inequality from stage n — 1, 

^n-l\ 1 



n-l 

E(-)^ 

fe=0 



to get 



k / qn~ 



E( 

fc=0 



-fc-1 



E(- 

fc=0 



\k-l 



'n-l\ 1 



k-l q, 



> 



n—k 



'n-l\ 1 



> 0. 



Qn-k 



This is the inequality of stage n for m = n — 1. (We have used the identity Q) = ("^^) + 
("Zi)-) Adding to it the inequality of stage n — 1 with m = n — 2 yields the inequality of 
stage n for m = n — 2. Repeating this process will give all the inequalities of stage n. □ 
It is helpful to introduce the quantities 

' (3.11) 
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Obviously, we have to = 1 (since qo = 1), and we have seen that the full set of inequalities 
restricting the Qn will be satisfied iff t.„ > for all n. (Recall we are assuming g„ > 0, Vn.) 



Moreover, given the tn, we can recover the Qn by inverting (3.11 



Lemma 4 




(3.12) 



Proof: This follows immediately from the identity 



k=0 



n—k I 



which itself follows from Equations (27) and (28) on the top of pg. 37 of 1 60 1. □ 
Thus, the t„ may be treated as free parameters (subject only to t„ > and to = 1)) a-^d the 
Qn can then be derived from ( |3.12| ). If this is done, the remaining transition probabilities 
an can be re-expressed more simply in terms of the t„, by inserting ( |3.12 ) into ( p.8| ) to get 



7^ = E*^E(-)'f 

'i'^ I k \ 



m\ (w — k 



w — m 
. w — I I 



whence 



l=in \ zu-l 



(3.13) 



Here, we have used an identity for binomial coefficients that can be found on page 63 of 
[27|. To express (|3.13| ) more concisely, define 



\(w,m) = ^ 



l=m. 




Then qn = A(n, 0) ^ and 



Or 



A(t37, m) 
A(n,0) 



(3.14) 



In this way, we arrive at the general solution of our inequalities. (Actually, we go slightly 
beyond our "genericity" assumption that an 7^ if we allow some of the tn to vanish; but 
no harm is done thereby.) 
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Let us conclude this section by noting that ( |3.12| ) implies 



g'o = 1 > gi > g'2 > g'3 > • • • (3.15) 

If we think of the g„ as the basic parameters or "coupling constants" of our sequential 
growth dynamics, then it is as if the universe had a free choice of one parameter at each 
stage of the process. We thus get an "evolving dynamical law", but the evolution is not 
absolutely free, since the allowable values of g„ at every stage are limited by the choices 
already made. On the other hand, if we think of the tn as the basic parameters, then the free 
choice is unencumbered at each stage. However, unlike the the t„ cannot be identified 
with any dynamical transition probability. Rather, they can be realized as ratios of two 
such probabilities, namely as the ratio Xn/qn, where Xn is the transition probability from 
an antichain of n elements to the timid child of that antichain. (Thus, if we suppose that 
the evolving causet at the beginning of stage n is an antichain, then tn is the probability 
that the next element will be born to the future of every element, divided by the probability 
that the next element will be born to the future of no element.) 



3.3.4 Proof that this dynamics obeys the physical requirements 

To complete our derivation, we must show that the sequential growth dynamics given by 
ni) or ([3.14]) obeys the four conditions set out in Section 3.2. 



3.3.4.1 Internal temporality 



This condition is built into our definition of the growth process. 
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3.3.4.2 Discrete general covariance 

We have to show that the product of the transition probabihties associated with a 
labeling of a fixed finite causet C is independent of the labeling. But this follows immediately 



from ( |3.8D [or ( 3.14 )] once we notice that what remains after the overall product 



|chi 

n 

j=0 

is factored out, is a product over all elements x £ C of poset invariants depending only on 
the structure of past(x). 

3.3.4.3 Bell causality 

Bell causality states that the ratio of the transition probabilities for two siblings depends 



only on the union of their precursors. Looking at (3^), consider the ratio of two such 
probabilities a„i and a„2- The g„ factors will cancel, leading to an expression which depends 
only on wi, W2, mi, and m2. Since these are all determined by the structure of the precursor 
sets, Bell causality is satisfied. 

3.3.4.4 Markov sum rule 

The sum rule states that sum of all transition probabilities an from a given parent C (of 
cardinality |C| = n) is unity. Since a child can be identified with a partial stem of the 
parent, we can write this condition, in view of ( p.l4| ), as 

i:i:.('f_-:3^i:.C) 

where S ranges over the partial stems of C. This must hold for any ti, since they may be 
chosen freely. Reordering the sums and equating like terms yields 



s 
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an infinite set of identities which must hold if the sum rule is to be satisfied by our dynamics. 
The simplest way to see that ( p^ is true is to resort to transitive percolation, for 



which ti = where t = p/q = p/{l — p) (c.f. (3.18) below). In that case we know that the 



sum rule is satisfied, so by inspection of ( |3.16| ), we see that the identity ( 3.17| ) must be true 



A more intuitive proof is illustrated well by the case of Z = 3. Group the terms on the 
left side according to the number of maximal elements: 

E (1^0°) + E m) + E + E (ifa = (?) 

S\m{S)=0 S\m{S)=l S\m{S)=2 S\m{S)=3 

+ E ('V) + E (1^1-2) + E 1 = Q 

S\m{S)=l S\m{S)=2 S\m{S)=3 

The first term is zero because the only partial stem with zero maximal elements is empty (i.e. 
\S\ =0). The second term is a sum over all partial stems with one maximal element. This 
is equivalent to a sum over elements, with the element's inclusive past forming the partial 
stem. The summand chooses every possible pair of elements to the past of the maximal 
element. Thus the second term overall counts the 3-element subcausets of C with a single 
maximal element. There are two possibilities here, the three-chain jand the "lambda" A. 
The third term sums over partial stems with two maximal elements, which is equivalent 
to summing over 2 element antichains, the inclusive past of the antichain being the partial 
stem. The summand then counts the number of elements to the past of the two maximal 
ones. Thus the third term overall counts the number of three element subcausets with 
precisely two maximal elements. Again there are two possibilities, the "V" V , and the 
. Finally, the fourth term is a sum over partial stems with three maximal elements, 
and this can be interpreted as a sum over all three element antichains • • • . As this example 
illustrates, then, the left hand side of ( p.l7| ) counts the number of I element subcausets of 
C, placing them into "bins" according to the number of maximal elements of the subcauset. 
Adding together the bin sizes yields the total number of / element subsets of C, which of 
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course equals C^). 

3.3.5 Sample cosmologies 

The physical consequences of differing choices of the t„ remain to be explored. To get an 
initial feel for this question, we list some simple examples. (Recall our convention that 
io = 1, or equivalently, qo = 1, where qo is the probability that the universe is born at all.) 

• "Dust universe" 

to = l, ti = 0, i>l 

This universe is simply an antichain, since, according to ( ^.12 ), qn = 1 for all n. An 



"example" is shown in Figure 3.8 



Figure 3.8: 16 element dust universe 



"Forest universe" 



tQ = ti = l; ti = 0, i>2 



This yields a universe consisting wholly of trees, since (see the next example) t2 = 
^3 = ^4 = ■ ■ ■ = implies that no element of the causet can have more than one past 
link. The particular choice of ii = 1 has in addition the remarkable property that, as 



follows easily from ( 3.14 ) , every allowed transition of stage n has the same probability 



l/(n + 1). Figure 3.£ displays an example. 




Figure 3.9: 23 element forest universe 
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Case of limited number of past links 



tj = 0, i > no 



Referring to expression ( 3.14| ) one sees at once that an vanishes if m > no- Hence, 
no element can be born with more than hq past links or "parents". This means in 
particular that any realistic choice of parameters will have t„ > for all n, since 
an element of a causal set faithfully embeddable in Minkowski space would have an 
infinite number of past links. See Figure 3.10| for an example. 



Figure 3.10: 9 element universe with ti = t2 = 1 



Transitive percolation 



tn. — t 



(3.18) 



We have seen that for transitive percolation, g„ = g", where q = 1 — p. Using 




Figure 3.11: 20 element transitive percolation universe with p = 1/4 
the binomial theorem, it is easy to learn from ( p. 12 ) or (3.11) that this choice of % 
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corresponds to tn = t"" with t = p/q. Clearly, t runs from to oo as p runs from to 



1. Figure 3.11 shows an example. 



A more lifelike choice? 

in = -7 (3.19) 
n! 

Due to its homogeneity (c.f. § |2.2.2| ), we have seen that transitive percolation with 
constant p yields causets which could reproduce — at best — only limited portions 
of Minkowski, de Sitter, or anti-de SitterP) space. It also suffers from a sort of "scale 
dependent dimension" which is incompatible with any continuum spacetime. This 
occurs because, at finite p, transitive percolation generates causets with approximately 
"constant width" (~ 1/p), meaning that the Hasse diagram looks roughly like an 
infinite cylinder with cross sectional area of roughly 1/p. Thus, at larger and larger 
length scales, the transitive percolation universe looks more and more like a one 
dimensional spacetime. Both of these issues suggest that one would have to scale p 
so that it decreased with increasing n. This implies that t„ should fall off faster than 
in any percolation model, hence (by the last example) faster than exponentially in 

2 

n. Obviously, there are many possibilities of this sort (e.g. t„ ~ e~"" ), but one 
of the simplest is tn ~ c/n!. This choice for the t„ would be our candidate of the 
moment for a most physically realistic choice of parameters. (The factor of t" in the 



numerator, which does not affect the asymptotic behavior, was added by Dou |26|, 
to simplify the study of the cosmological renormalization behavior of this dynamics. 
In fact, under the renormalization fiow, this dynamics approaches that of originary 
transitive percolation with a t which approaches 0. A 20 element sample causal set 



^^We have already seen that transitive percolation cannot yield all of any homogeneous space, due to the 
presence of "posts". It will fail to produce all of anti-de Sitter for yet another reason — that intervals of 
sufficiently large size do not have finite volume, so that no locally finite order can embed faithfully into such 
a spacetime. 
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"generated by" this dynamics is shown in Figure 3.3. 



3.4 Originary dynamics 

It is possible that we can arrive at every possible solution of our conditions by taking limits 
of the generic solution. One such special theory is the originary percolation model, which 



was introduced in §2.2.5. It is the same as the transitive percolation model, but with the 
added restriction that each element except the original one must have at least one ancestor 
among the previous elements. Algorithmically, we generate potential elements one by one, 
exactly as for plain percolation (by "joining" to each existing element with probability p, 
then adding relations required by transitivity) but discard any such element which would 
be unrelated to all previous elements. Causets formed with this algorithm always have a 
single minimal element, an "origin". The transition probabilities for originary percolation 
are just those of ordinary transitive percolation with an added factor of (1 — g") in the 
denominator at stage n. 

This can be generalized to a non-percolation dynamics. Here the causal set grows as 
usual for a general dynamics, but with the added restriction that it must always possess a 
1-chain as a full stem. We call this an originary dynamics. The poset of originary causets 
can be transformed into the poset of all causets (exactly) by removing the origin from every 
originary causet. 

Further generalizations are also possible, in which a more complex stem of the causet is 
enforced, e.g. the restriction that after the first two elements form, the causet must always 
have a 2-chain as a full stem (or partial stem). However, not every poset can be used as a 
stem in this manner, for many choices are incompatible with Bell causality. 

We conjecture that each of these exceptional families of solutions are singular limits of 
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the generic family. For example, originary percolation is the A — > oo limit of the dynamics 
given by to = 1, tn = Af^, n = 1, 2, 3, . . . 



3.5 The stochastic growth process as such 

We have seen that, associated with every labeled causet C of size N, is a net "probability of 
formation" Pr(C) which is the product of the transition probabilities Oi of the individual 
births described by the labeling: 

Af-l 

Pr(C) = l[ai (3.20) 



i=Q 



where ai = a{i,Wi,mi) is given by ( |3.8| ) or ( |3.14| ), Wi and rrii are respectively the size and 
number of maximal elements in the precursor (= past) of the element labeled i. Using 
( 3.14| ), we can write this more explicitly as 



nf=o'A(i,0) 



Pr(C) = ^^'=°„;\7'::^ . (3.21) 



We have also seen that Pr(C) is in fact independent of the labeling, i.e. Pr(C) = Pr(C') 
where C is the same causet as (7, but with a different labeling. 

The net probability of arriving at an unlabeled causet C, at stage of the growth 
process is 

Pr7v(C) = W{C) Pr(C) (3.22) 

where C is the causet C endowed with any (natural) labeling, N = |C|, and W{C) is 
the number of inequivalentj^ labelings of C, or in other words, the total number of paths 
through V that arrive at C, each link being taken with its proper multiplicity. Expressing 



^^Two labelings of C are equivalent iff they are related by an automorphism of C. 
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(3.21) more intrinsically, we can write (p.22|) as 



n X{'uj{x),m{x)) 

Fin{C) = WiC) ^^^^-^ , (3.23) 

n A(i,0) 

j=0 

where tu(x) = |past2;| and m{x) = |maximal(past2;)|. This expression, as far as it goes, is 
manifestly "covariant" in the sense explained above. Also causality is manifest in the sense 
that it is expressed as a product of factors, one for each element (save the W{C)), each of 
which depends only on that element's past. However, as explained in §3.2.2 and § |1.3.1| , it 



has no direct physical meaning. Here we briefly discuss some probabilities which do have a 
fully covariant meaning and show how, in simple cases, they are related to ^ oo limits 
of probabilities like ( 3.23| ) . 



As a rudimentary example of a truly covariant question, let us take "Does the two-chain 
ever occur as a partial stem of C?". The answer to this question will be a probability, P, 
which it is natural to identify as 

P= lim PrAr(X7v), (3.24) 

where Xj^ is the event that "at stage A^" , C possesses a partial stem which is a two-chain. To 
state this more precisely, define to be the set of A^-orders which satisfy some criterion, 
e.g. that the order possess the causet as a partial stem (5* was a two-chain in the above 
example), and define 

Ptn{Xn)= E Prjv(C). (3.25) 

Of course it is not guaranteed that P will be well defined in the N ^ oo limit. In this 
connection, we conjecture that the questions "Does S occur as a partial stem of C?" furnish 
a physically complete set, when S ranges over all (isomorphism equivalence classes of) finite 
causets. 
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As a simple example of an answer to a question of this form, consider "What is the 
probability P that a 1-chain is a full stem of the universe C?". Clearly this is equivalent 
to the question "Does C have a unique minimal element?". In terms of partial stems, this 
question is equivalent to demanding that the 2-antichain not be a partial stem of C. The 
answer is simple to formulate by thinking in terms of the growth process, as follows. At 
stage of the process, C is a 1-chain, which occurs with probability 1. At stage 1, C must 
not become a 2-antichain, which occurs with probability l — qi- At stage 2, the new element 
must not be born with no ancestors, which occurs with probability 1 — ^2- At stage 3, the 
same condition occurs with probability 1 — ^3, and so on. Thus in the limit N ^ oo, the 
answer becomes 

oo 

= 11(1 - Qi) . (3.26) 

1=1 

Expressed in terms of the tn, the qi in ( p.26| ) are simply replaced with l/J2k=o {k)^k- 



3.6 Two Ising-like state-models 

The dynamics, as written in ( |3.20| ) (say) can be expressed in terms of either of a pair of 



Ising-like state models. To derive one such model, consider a„, as given in (3^). Inserting 
this form into ( 3.20D (and discarding the labeling decoration on C) 



N-l rrij / \ 1 

j=0 k=o V / 

where we have placed an index j on zu and m to indicate that they refer to the transition 
at stage j. For a given j, the sum J^t^o i"k) regarded as a sum over subsets of the 

nij maximal elements to the past of element j (where k can be regarded as the cardinality 
of each subset), or equivalently a sum over the nij links whose future endpoint is j. This 
suggests an interpretation of this as a sum over Z2 valued "spin configurations" on these 
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links. Before describing in detail how this construction follows, it will be helpful to state a 
number of definitions. 

The following paragraph will refer to a given causal set C. First define TZ as the set of 
all relations in C 

n = {{x,y) C X C\x <y} . 

TZj will denote the set of relations in C whose future endpoint is the element labeled j 

TZj = {{x, Xj) \ X ^ Xj y X £ C} . 

Now define 4> as a. map which assigns a Z2 (= {0, 1}) valued "spin" to each relation of the 
causal set C, i.e. 

: 7^ ^ Z2 . 

Let the set of all such maps, for the given causal set C, be denoted by Furthermore, 
define a restriction (pj of a map (p to an element j by restricting the domain to only those 
relations which have j as a future endpoint, i.e. 

We denote the set of such restricted maps by ^j. We also need to define two quantities 
associated with a map (f>, a{cj)) and r{cj)). The former counts the number of "absent spins" 
in (/>, i.e. it represents the number of relations which map to 

a((/>) = \{x G 7^|0(x) = 0}| . 

The latter counts the number of "present spins" in (p, i.e. 

r[(P) = \{x G 7^|(/)(x) = 1}| . 

It is to be understood that a restricted map (pj can be used in the place of with TZ 
replaced by TZj in the above two definitions. 
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In order to express ( |3.27| ) using these definitions, it is necessary to place a restriction on 
all maps cf) introduced in the previous paragraph, namely that each relation which is not a 
link (i.e. pairs {x,y) such that int[x,y] 7^ 0) map to 1. Such maps (and sets of such maps) 
will be decorated with a With these definitions in place, ( |3.27| ) can be written as 

i=o / i=o r.pj; %{^,) 



where we have interpreted the index k in ( 3.27| ) as counting the number of "zero spins" on 



links "pointing to" j. In writing (3.28), we have noted that Wj, the number of elements to 



the past of j, counts the number of relations in the domain of Thus, after subtracting 
the "zero spins", the subscript in which it appears becomes simply r{(j)j). 

Equation ( |3.28| ), save an initial coefficient (jjfSQ^Qj^, is written as a product, one 
factor for each element j, of sums of terms, one for each "spin configuration" at j. (A "spin 
configuration at j" being an assignment of ones or zeros to each past-directed link at j.) 
Expanding, we arrive at a sum of terms, each of which contains one factor from each element 
j, for one choice of spin configuration at j. The sum contains a term for each possible spin 
configuration at each j. Thus, after expanding, (|3.28| ) becomes 

/n~i \ N-i _ 1 



or, using a more covariant notation. 



'N-l 



p^ic) = n E n i-r^""-^-^ ■ (3.29) 



Equation ( 3.29| ) writes the probability of arriving at a particular causal set C, at stage 



N of the growth process, as a sum over spin configurations on the causal set, where only the 
spins on the links are permitted to vary. For each such spin configuration, each element x of 
C contributes a "vertex factor" of (— )"(''^^) ^ . If these vertex factors are to be interpreted 
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as Boltzman weights, then the negative values for odd numbers of "present" past-hnks are 
a bit pecuhar. 

A second model arises by inserting ( p. 13 ) into (|3.20| ). 



N-l y-^i f^j-mjU, 

p.-(c) = n t:] X (3.30) 



From ( p. 12 ) the factors in the denominator are easily seen to be simply the overall product 



Y{^=Q l/^j)- As before, the sum in the numerator, YHilm - can be interpreted 

as a sum over subsets of relations. This time, however, the sum is over subsets of relations 
which are not links. To express this in terms of spin configurations 4>, constrain all such 
maps to yield 1 on links, but allow them to vary freely on the non-link relations. We 
decorate maps which respect such a constraint with a Then 

n n E Kify > 

where we have interpreted the index j in ( |3.30| ) as the number of "present" relations in 
each 4>j. Expanding this product, as before, leads to a sum of terms, one for each "spin 
configuration" on the entire causal set C, of a product of factors, one for each element. 
Thus 



'N-l \ N-l 



or, in a more covariant notation 



'N-l 



^^(c) = I n E n Ki,Z) ■ (3-31) 

J=0 1 Jg$a;eC 



Equation (3.31 ) represents a second way to express the probability of arriving at a causet 



C in terms of a model of Ising-like spins on its relations. Here the spins on the links are 
fixed at 1, while those on the other relations are free to vary. This time all vertex factors are 
positive, in closer agreement with what one would expect from physical Boltzmann weights. 
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These two models (and especially the second) show that the sequential growth dynamics 
can be viewed as a form of "induced gravity" obtained by summing over ("integrating out") 
the values of underlying spin variables a. This underlying "matter" theory may or may not 
be physically reasonable (Does it obey its own version of Bell causality, for example?), 
but at a minimum, it serves to illustrate how a theory of non-gravitational matter can be 
hidden within a theory that one might think to be limited to gravity alone.|^ |^ It should 
be noted that these "spin models" are "non-interacting" in that each "lattice site" has its 
own "reserved" set of spins which affect the value of only its vertex factor, with no two 
lattice sites "sharing" any spins. In order for these spin models to give non-trivial results, 
an effective interaction must emerge from the gravitational dynamics in the sum over causal 
sets. 

3.7 Further Work 

The sequential growth dynamics can be simulated directly on a computer, but only for very 
small N. For t„ = 1/n! it takes a minute or so to generate a 64 element causet on a DEC 
Alpha 600 workstation. Because the number of partial stems, and hence the number of 
possible precursors for a new element, of a element causal set grows like 2^, it is difficult 
to simulate the growth process directly. A workaround may involve using something like a 
metropolis algorithm at each stage to select a precursor. Issues such as detailed balance for 

^•^In this connection, it bears remembering that Ising matter can produce fermionic as well as bosonic 
fields, at least in certain circumstances. |46| 

^*References and describe a similar example of "hidden" matter fields in the context of 2- 
dimensional random surfaces (Euclidean signature quantum gravity) and the associated matrix models in 
the continuum limit. Unfortunately, the matter fields used (Ising spins or "hard dimers") were unphysical 
in the sense that the partition function was a sum of Boltzmann weights which were not in general real and 
positive. This is much like our first state model described above. To the extent that the analogy between 
these two, rather different, situations holds good, our results here suggest that there might be, in addition 
to the matter fields employed in |^^, another set of fields with physical choices of the coupling constants, 
which could reproduce the same effective dynamics for the random surface. 
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stepping through precursors would have to be sorted out. 

Analytic results, so far, are available only for the special case of transitive percolation. 
An important question, of course, is whether some choice of the tn can reproduce general 
relativity, or at least reproduce a Lorentzian manifold for some range of t's and of n = \C\. 



Another question is whether the "Ising matter" introduced in §3.(: gives rise to an 



interesting effective field theory, and what relation it has with the local scalar matter on a 



background causal set studied in 52] 



Another possibility for obtaining the behavior of a scalar field on a causal set arises 
if we are willing to drop the acyclicity property of the order (i.e. replace the transitive. 



irreflexive order with a transitive, reflexive "pseudo-order"), as mentioned in §1.2.5. This 
relaxation allows the possibility for causal cycles to exist, but with the property that each 
element of a cycle has the same causal relations with the remainder of the pseudo-order 
as every other point in the cycle. Thus at the level of the pseudo-order, these cycles hold 
no more information than that of one element in an ordinary partial order, along with a 
positive integer representing the "degeneracy" of the cycle. This indicates that a pseudo- 
order is equivalent to an order with a positive integer attached to each element, where each 
integer represents the size of a cycle which exists at that point. A generalization of our 
dynamics to allow such a possibility may be achieved simply by allowing, for a causal set 
with m maximal elements at some stage of the growth process, m additional transitions, 
corresponding to incrementing the integer attached to one of the maximal elements by 1. 

Another set of questions concerns the possibility of a more "manifestly covariant" formu- 
lation of our sequential growth dynamics - or of more general forms of causal set dynamics. 
Can Bell causality be formulated in a gauge invariant manner, without reference to a choice 
of birth sequence? Is our conjecture correct that all meaningful assertions are logical com- 
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binations of assertions about the occurrence of partial stems/past sets? 

Also, there are the special cases we left unstudied, for example the originary dynamics 
and its generalizations. Are there other special, non-generic cases of interest? In addition, it 
may be interesting to see what sort of dynamics may arise from omitting the combinatorial 
coefficients in the Markov sum rule, which corresponds to imparting an indistinguishability 
condition on the individual causal set elements. Such a condition is perhaps not unreason- 
able physically. 



Chapter 4 



Conclusions 

4.1 Summary 

The Causal Set approach to quantum gravity selects a very sparse framework for the "sub- 
stance" of the theory, and seeks to "recover" most of existing physical theory in some ap- 
propriate limit. A disadvantage of this philosophy is the difficulty it raises by abandoning 
many of the "typical" structures and methods that are familiar to us. In a way this works 
itself out as an advantage, however, in that it forces us to think critically about many issues 
at the foundation of physical theory. In doing so, the structure which arises is robust in the 
sense that issues such as Lorentz invariance, background independence, general covariance, 
discreteness, non-perturbative formulation, Lorentzian signature, have been addressed at 
the outset, rather than being put off until the theory is more developed. Also the clarity of 
the fundamental structure provides a fertile ground for addressing philosophical issues such 
as the nature of quantum causality in a closed system. 

Another unique aspect of Causal Sets is the theory's departure from some conventional 
intuition in a number of respects. Most obvious is the rejection of the continuum as funda- 
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mental. Another is the complete departure from determinism — even at the fundamental 
scale, when all aspects of the "state of the system" are "known" , the classical limit of the 
theory is postulated to be stochastic in nature. (As opposed to the philosophical attitude 
of kinetic theory, which assumes that only the incomplete knowledge of the state leads to 
indeterminism, c.f. [^4|) Unitarity will likely have to be abandoned to formulate the quan- 
tum theory in a discrete setting. Qualitatively, the theory predicts a non-zero cosmological 



constant |57, 43|. Locality as a fundamental physical principle seems to be abandoned. 
The problem here essentially extends from the fact that the links, which are the analogs of 
nearest neighbors in a Euclidean signature lattice, extend arbitrarily far into the past and 
to arbitrary spatially distant points in some reference frame. The work of |23, 52 1 indicates 
how effective locality is preserved in the causal set, but, when looking at the discrete order 
alone, it arises in a complicated and non-intuitive manner. 

Before beginning this work, the Causal Set theory was at a stage where a lot was known 
about kinematical issues, but there seemed to be many obstacles to the construction of a 
dynamics for the theory. Even if some reasonable guess for the action was made, the issue 
of how to do the sum over histories to compute the measure was difficult, because of the 
sum over an enormous discrete sample space. This led to the search for an algorithm to 
sample the space of causal sets. The proposed algorithm, transitive percolation, failed to 
perform the desired task, but did suggest itself as a toy dynamics. As a dynamical model, 
it has a number of promising features, and the possibility that it reproduces a region of 
continuum spacetime has not yet been excluded. 

By thinking in terms of a stochastic growth process, and positing some very basic 
principles, we were led almost uniquely to a family of dynamical laws (stochastic processes) 
parameterized by a countable sequence of coupling constants g„, (or equivalently the 
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This result is quite encouraging in that we now know how to speak of dynamics for a theory 
with discrete time. In addition, these results are encouraging in that there exists a natural 
way to transform this theory, which is expressed in terms of a classical probability measure, 
to a quantum measure. A sketch as to how one might proceed in doing this is provided 
below Thus there is good reason to expect that we are close to constructing a background 
independent theory of quantum gravity. 



4.2 Quantum Dynamics 

Since our theory is formulated as a type of Markov process, and since a Markov process 
mathematically is a probability measure on a suitable sample space, the natural quantum 
generalization would seem to be a quantum measure (or equivalent "decoherence func- 
tional") on the same sample space. The question then would be whether one could find 
appropriate quantum analogs of Bell causality and general covariance formulated in terms 
of such a quantum measure. If so, we could hope that, just as in the classical theory treated 
herein, these two principles would lead us to a relatively unique quantum causal set dy- 
namicSjQ or rather to a family of them among which a potential quantum theory of gravity 
would be recognizable. Let me sketch briefly how one might go about constructing this 
quantal generalization. 

The quantum dynamics for causal sets will be expressed in terms of a quantum measure, 
which is a generalization of a classical probability measure. For a more detailed discussion, 
see 1 58, 56, 51]. It is helpful to express the measure in terms of a decoherence functional 



D{C' , C") which assigns a complex number to pairs of histories C and C" . In the context 
of causal sets, a "completed" causal set C is regarded as a history. (Recall that a completed 



^See for a promising first step toward such a dynamics. 
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causal set is one which has infinite cardinahty; it has "run to completion".) More correctly, 
the decoherence functional will be defined for pairs of sets of histories, which has the 
following properties: (for any (disjoint) sets of histories R, S, T) 

positivity: D{S,S) > 

additivity in each argument: D{R U S,T) = D{R, T) + D{S, T) (4.1) 
hermitian: D{R, S) = D{S, R)* (4.2) 

(U indicates the union of disjoint sets.) The quantum measure of a set S, fi{S) is given by 
the diagonal elements of D, i.e. //(S*) = D{S,S). Two additional properties are that for a 
set S of measure zero fi{S) = 0, D{A, 5) = for any set A, and a quantum generalization 

of (iHD 

H{AUBUC)- fi{A UB)- fi{A UC)- fi{B UC) + n{A) + fi{B) + fi{C) = . 

In practice, since the expressions we know how to write down for causal set "amplitudes" 
are written in terms of finite causal sets, as mentioned earlier, we conjecture that the 
appropriate sets to consider are cylinder sets, defined by specifying a labeled causet as a 
partial stem, and including in that set all possible extensions of that labeled stem into the 
future. Thus with each pair of labeled finite causal sets C' , C" we can associate a complex 



number D{C', C"), which may look like some quantum generalization of ( 3.21 ). 

To construct a dynamics for such an object, this language needs to be re-expressed in 
terms of "transition amplitudes" , so that the growth process construction can be carried over 
to the quantum dynamics. This may be accomplished by defining a transition amplitude 
T{C", C C) by 

rr(r<" n r<i\ _ ^i^" ,C') 
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(Note that there is some danger here if we aUow some of the D to vanish, which should not 
be uncommon in the quantum theory due to complete destructive interference, but we'h 
evade this issue for now. It is hkely that this problem can be overcome using limits.) 

With this definition in place, to derive the quantum dynamics, we need to generalize 
our physical conditions of ^3.2| . Internal temporality should carry over directly, since it 
is encoded into the definition of the growth process itself. General covariance retains its 
obvious meaning, that products of transition amplitudes leading to a given finite causal set 
in the sequential growth process be independent of labeling. (This translates into the state- 
ment that D{A,B) is invariant under relabeling of its arguments.) The Markov sum rule 
generalizes to additivity in each argument (4.1). For Bell causality, a promising generaliza- 



tion is to demand that for one argument S fixed, D(S,-) obeys the classical Bell causality 
condition. This version of quantum Bell causality seems to be obeyed by non-relativistic 
quantum mechanics, assuming that there are no correlations existing in the initial data. 

Phrased in this manner, this definition for quantum causality makes the quantum theory 
look a lot like "the classical measure, squared". However, the hermitian condition ( [4.2[ ) 
provides a constraint on the theory which may considerably limit this freedom. The hope is 
that one can follow the general methodology used in deriving Pr(C) to construct a quantum 
dynamical law for causal sets, phrased in terms of the D{C', C"). 



Appendix A 

Consistency of physical conditions 



Our analysis of the conditions of Bell causality et al. unfolded in the form of several lemmas. 
Here we present some similar lemmas which strictly speaking are not needed in the present 
context, but which further elucidate the relationships among our conditions. We expect 
these lemmas can be useful in any attempt to formulate generalizations of our scheme, in 
particular quantal generalizations. 

Lemma 5 The Bell causality equations are mutually consistent. 



Proof: The top of Figure shows three children of an arbitrary causal set C„. The 
shaded ellipses represent portions of Cn- The small square indicates the new element whose 
birth transforms C„ into a causal set Cn+i of the next stage. The smaller ellipse "stacked on 
top of" the larger ellipse represents a subcauset of Cn which does not intersect the precursor 
set of any of the transitions being considered (i.e. none of its elements lie to the past of 
any of the new elements). This small ellipse thus consists entirely of "spectators" to the 



transitions under consideration. The bottom part of Figure A.l shows the corresponding 
parent and children when these spectators are removed. 
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union of precursors 



Figure A.l: Two families related by Bell causality 

Notice that one of the three children is the gregarious child. We will show that the 
Bell causality equations between this child and each of the others imply all remaining Bell 
causality equations within this family. Since no Bell causality equation reaches outside 
a single family (and since, within a family, the Bell causality equations that involve the 
gregarious child obviously always possess a solution — in fact they determine all ratios of 
transition probabilities except for that to the timid child), this will prove the lemma. 

In the figure ri and r2 represent a general pair of transitions related by a Bell causality 
equation, namely 



n _ £i 

r2 S2 



(A.l) 



But, as illustrated, each of these is also related by a Bell causality equation to the gregarious 
child, to wit: 



— - — and — - — 



n S3 



^3 S3 



(A.2) 
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Since (A.l) follows immediately from (A. 2), no inconsistencies can arise at stage n, and the 
lemma follows by induction on n. □ 

Lemma 6 Given Bell causality and the further consequences of general covariance that are 
embodied in Lemma |^ all the remaining general covariance equations reduce to identities, 
i.e. they place no further restriction on the parameters of the theory. 



Proof: 




Figure A. 2: Consistency of remaining general covariance conditions 



Discrete general covariance states that the probability of forming a causet is independent 
of the order in which the elements arise, i.e. it is independent of the corresponding path 
through the poset of finite causets. Now, general covariance relations always can be taken 
to come from 'diamonds' in the poset of causets, for the following reason. As illustrated in 



Figure A. 2, any pair of children of a causet (siblings) will have a common child obtained by 
adjoining both new elements of the two siblings, i.e. adding to the "grandparent" both the 
new element which defines one sibling and the new element which defines the other sibling. 
(For example, consider the case where the 2-antichain is the grandparent and it has the 
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child ••• (by adding a disconnected element) and the child (by adding an element to the 
future of both elements of • • ) . To find their common child A. add a disconnected element 
to A, or an element to the future of two of the elements of ••••) Now, still referring to 



Figure A. 2, let \D\ = n and suppose inductively that all the general covariance relations 
are satisfied up through stage n. A new condition arising at stage n -\- 1 says that some 
path arriving at B via x has the same probability as some other path arriving via z. But, 
by our inductive assumption, each of these paths can be modified to go through D without 
affecting its probability. Thus, the equality of our two path probabilities reduces simply to 
ax = bz. 

Now by Bell causality and lemma 

X b 



whence 



qn q-n-i 



ax = ab- 



qn-i 



But by symmetry, we also have 



therefore ax = bz, as required. □ 



bz = ba ^" 



In-l 
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